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ABSTRACT 


Of  the  physical  effects  that  limit  the  propagation  of  a high  intensity 
laser  beam  through  the  atmosphere,  air  breakdown  has  the  most  catas- 
trophic influence.  Completely  ionized  air  attenuates  10.  6 ^ radiation  within 
a few  centimeters.  It  is  therefore  imperative  that  the  threshold  for  air 
breakdown  under  realistic  circumstances  be  determined.  In  this  report 
the  physics  of  long  pulse  laser-induced  breakdown  process  in  particle- 
contaminated  air  is  studied  in  four  stages.  The  first  stage  is  particulate 
heating  and  vaporization.  Simple  models  are  used  to  predict  the  time  to 
vaporization  at  the  particle  surface,  the  time  uo  vaporize  a given  fraction 
of  the  particle,  and  the  pressure  of  the  vapor  immediately  above  the 
particle  surface.  The  second  stage  is  the  heating  of  the  vapor  cloud.  The 
heating  rate  is  calculated  under  two  approximations.  For  a short  time 
period  after  the  vapor  cloud  is  formed,  its  central  region  is  assumed  to  be 
effectively  insulated  from  the  air,  and  thermal  conduction  losses  can  be 
neglected.  At  later  times  thermal  conduction  losses  are  important  through- 
out the  whole  ..  oud;  then  a model  in  which  the  conduction  losses  are 
averaged  over  the  volume  of  the  vapor  is  used.  Using  these  models  we 
are  able  to  estimate  the  time  for  the  vapor  cloud  to  reach  a high 
temperature  (nominally  chosen  to  be  20,000°K)  and  the  minimum  laser 
intensity  required  to  heat  the  vapor  cloud.  The  third  stage  of  the  break- 
down process  is  the  creation  of  air  plasma.  An  estimate  is  made  of  the 
tamperature  the  surrounding  air  must  reach  in  order  to  readily  absorb 
the  laser  energy.  Then,  from  crude  estimates  of  the  energy  flux  into  the 
air,  the  rate  of  expansion  of  the  plasma  into  the  air  is  determined.  The 
final  stage  is  the  growth  of  the  air  plasma.  The  growth  is  followed  as  a 
function  of  time  for  various  values  of  the  initial  air  plasma  size  and  the 
laser  intensity. 
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1.  INTRODUCTION 


Of  the  physical  effects  that  limit  the  propagation  of  a high  intensity 
laser  beam  through  the  atmosphere,  air  breakdown  has  the  most  catastrophic 
influence.  Completely  ionized  air  attenuates  10.  6pi  radiation  within  a few 
centimeters.  It  is  therefore  imperative  that  the  threshold  for  air  break- 
down under  realistic  circumstances  be  determined. 

Studies  of  the  breakdown  threshold  in  clean  ai  r and  the  noble  gases 
are  numerous,  and  observed  thresholds  are  in  good  agreement  with  an 
extension  of  microwave  breakdown  theory  to  optical  frequencies.  When  air 
is  contaminated  with  particles,  it  has  been  observed  that  the  breakdown 
threshold  is  lowered  substantially  below  the  clean  air  value.  The  mecha- 
nisms by  which  particles  induce  breakdown  are  not  clearly  understood. 

However,  it  is  generally  accepted  that  breakdown  occurs  as  a result  of 
vaporization  of  the  particles,  heating  of  the  vapor  and  surrounding  air,  and 
ultimately  ionization  and  "breakdown"  of  the  air  within  the  laser  beam  volume. 
The  sequence  of  events  and  the  detailed  physical  processes  that  occur  is 
strongly  dependent  upon  laser  intensity,  wavelength  and  pulse  time. 

At  10.  6|im  and  for  short  laser  pulses  (<  1 Usee)  the  observed  breakdown 
threshold  is  one  to  two  orders  of  magnitude  lower  than  the  clean  air  value  for 
particle  sizes  from  3 |im  to  approximately  30  um  in  diameter.  The  threshold 
is  fairly  insensitive  to  the  particle  materials  that  have  been  studied  and 
decreases  with  increasing  particle  diameter.  In  addition,  for  a pulse  length  of 
approximately  100  nsec,  Lencioni2has  found  that  the  breakdown  threshold  in  the 
presence  of  particles  scales  as  the  inverse  square  of  the  laser  wavelength.  This 
is  consistent  with  cascade  breakdown  in  either  the  particle  vapor  or  the  air 

3 

surrounding  the  particle.  Triplett  and  Boni  have  performed  a theoretical 
study  of  the  interaction  of  laser  radiation  with  a single  suspended  particle. 
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Their  analysis  applies  only  when  the  laser  intensity  is  sufficiently  large  that 
the  rapidly  evaporating  or  " exploding"  particle  drives  a strong  shock  wave 
into  the  surrounding  air.  Breakdown  then  occurs  via  a nonequilibrium 
cascade  process  in  either  the  vapor  or  shock  heated  air.  Calculations 
indicate  that  reductions  in  the  breakdown  threshold  of  slightly  greater  than 
one  order  of  magnitude  below  that  of  clean  air  are  feasible  by  this  mechanism. 

An  experimental  study  of  laser-induced  breakdown  in  contaminated 

air  with  long  (100  Msec)  10.6  Urn  laser  pulses  was  performed  by  Schlier, 
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Pirri  and  Reilly.  The  most  significant  result  of  this  study  was  that  the 

laser  intensity  for  breakdown  in  contaminated  air  was  almost  three  orders 

of  magnitude  below  the  clean  air  value.  This  result  is  one  to  two  orders 

of  magnitude  below  experimental  and  theoretical  predictions  for  breakdown 

with  microsecond  or  submicrosecond  pulses.  In  addition,  an  "incubation" 

time  or  delay  time  between  laser  turn  on  and  breakdown  (as  dr  cted  by  a 

reduction  in  the  transmitted  laser  power)  was  observed.  Luminous  plasmas 

were  observed  during  the  incubation  period.  Often  a plasma  was  seen  but  no 

breakdown  occurred.  There  is  at  present  no  theoretical  understanding  of  the 

contaminated  air  breakdown  mechanisms  for  these  long  laser  pulses  (>  lOpsec). 

Since  the  observed  delay  times  are  10  - ?0psec,  the  breakdown  mechanism  cannot 

be  nonequilibrium  cascade  ionization  of  shock  heated  air  or  vapor.  It  appears 

that  the  vapor  formed  by  the  evaporating  particle  heats  in  local  thermodynamic 

equilibrium  and  transfers  its  energy  to  the  surrounding  air  by  conduction  and 

radiative  transport.  The  air  ultimately  heats  via  inverse  Bremsstrahlung 

absorption  until  the  formed  plasma  becomes  opaque  to  the  laser  radiation. 

This  plasma  formation  mechanism  is  similar  to  that  which  leads  to  the  ignition 
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of  laser- supported  combustion  waves  above  surfaces.  ' The  wavelength 
dependence  of  long  pulse  thresholds  is  unknown  at  present. 

In  this  report  the  physics  of  long  pulse  laser-induced  breakdown  process 
in  particle  - contaminated  air  is  studied  in  four  stages.  The  first  stage,  which 
is  discussed  in  Section  2,  is  particulate  heating  and  vaporization.  Simple 


modal*  are  used  to  predict  the  time  to  vaporisation  at  the  particle  surface, 
the  time  to  vaporise  a given  fraction  of  the  particle,  and  the  pressure  of 
the  vapor  immediately  above  the  particle  surface.  The  second  stage,  the 
heating  of  the  vapor  cloud,  is  analysed  in  Section  3.  The  heating  rate  is 
calculated  under  two  approximations.  For  a short  time  period  after  the 
vapor  cloud  is  formed,  its  central  region  is  assumed  to  be  effectively 
insulated  from  the  air,  and  thermal  conduction  losses  can  be  neglected. 

At  later  times  thermal  conduction  losses  are  important  throughout  the 
whole  cloud;  then  a model  in  which  the  conduction  losses  are  averaged 
over  the  volume  of  the  vapor  is  used.  Using  these  models  we  are  able  to 
estimate  the  time  for  the  vapor  cloud  to  reach  a high  temperature  (nominally 
chosen  to  be  20,000°K)  and  the  minimum  laser  intensity  required  to  heat  the 
vapor  cloud.  The  third  stage  of  the  breakdown  process,  the  creation  of 
air  plasma,  is  studied  in  Section  4.  An  estimate  is  made  of  the  temperature 
the  surrounding  air  must  reach  in  order  to  readily  absorb  the  laser  energy. 
Then,  from  crude  estimates  of  the  energy  flux  into  the  air,the  rate  of  expan- 
sion of  the  plasma  into  the  air  is  determined.  The  final  stage,  the  growth 
of  the  air  plasma,  is  analysed  in  Section  5.  The  growth  is  followed  as  a 
function  of  time  for  various  values  of  the  initial  air  plasma  size  and  the 
laser  intensity. 

A plot  of  the  experimentally  determined  threshold  intensity  versus 

time  to  breakdown  is  shown  in  Fig.  1.  I for  10.  6n  radiation.*  The  region 
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of  interest  in  this  report  is  the  intensity  range  from  3x10  W/cm  to 
6 2 

8x10  W/cm  . Spherical  particles  of  having  a diameter  of  30|im 

are  observed  to  cause  breakdown  within  20  to  lOO^sec  from  the  beginning 

of  the  pulse.  Glass  fibers  also  lead  to  similar  breakdown  times,  but  they 

are  not  analysed  in  this  report.  For  comparison  with  data,  numerical 

calculations  are  made  in  each  section  for  AUO,  particles  30  nm  in  diameter 
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exposed  to  10.  6^  radiation  of  intensity  between  3 x 10  W/cm  to  8 x 10  W/cm 


THRESHOLD  INTENSITY  (W/cm 


-9  •«  *7  -6  -5  -4 

10  10*  10 r 10*  10°  10*  I 

PULSE  TIME -TIME  TO  BREAKDOWN  (sec) 


Fig,  1.  1 Measurement*  of  Air  Breakdown  in  the 
Presence  of  Particulate  Matter 
(From  Ref,  1) 


Du*  to  th«  complicated  coupling  between  heating,  vapor  and  air  dynamics 

simple  acaling  laws  for  the  variation  of  breakdown  time  and  threshold  as 

a function  of  laser  wavelength  and  particulate  composition  can  not  be  given. 

However,  in  each  stage,  the  important  parameters  of  the  breakdown  process 

which  change  with  wavelength  and  particulate  composition  are  identified. 

In  Section  6,  this  report  is  concluded  with  the  calculation  of  the  entire 
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breakdown  sequence  for  a laser  intensity  at  5 x 10  W/cm  . The  conclusions 
of  each  of  the  four  stages  are  summarised  and  experiments  are  suggested. 


2.  PARTICULATE  1EATING  AND  VAPORIZATION 


A particulate  exposed  to  a laser  beam  absorbs  energy  and  thereby 
increases  its  temperature.  If  the  beam  is  intense  enough  the  particle  may 
reach  the  vaporization  temperature.  Many  interesting  phenomena  can 
occur  during  the  process;  for  example,  the  particle  may  move  due  to  the 
pressure  of  the  vapor,  or  the  particle  may  shatter  because  of  unequal 
heating.  However,  for  the  purpose  of  estimating  breakdown  thresholds 
the  important  quantities  are  the  size  and  composition  of  the  vapor  cloud 
formed  by  the  particle.  In  this  section,  therefore,  simple  models  are 
employed  to  characterize  (1)  the  time  when  bulk  vaporization  begins, 

(2)  the  time  required  to  vaporize  an  arbitrary  fraction  of  the  particle,  and 

(3)  the  pressure  of  the  vapor  above  the  particle  surface. 

The  optical  properties  of  the  particulate  are  determined  by  its 
index  of  refraction  which  has  both  a real  part  nj  and  an  imaginary  part 
n ; and  by  the  ratio  of  the  radiation  wavelength  \ to  the  radius  R 

M 

of  the  particulate.  For  R « X the  Rayleigh  limit  for  the  absorption  cross 

section  is  approached.  For  large  particles,  R » X,  geometric  optics  applies. 

If,  in  addition,  the  complex  index  of  refraction  n obeys  | n - 1 | « 1 and 

4TT  n /RX  < 1,  the  cross  section  is  given  by  volume  absorption  with  the  bulk 
2 

absorption  coefficient  a defined  by 

a = 4nn2/X  (2.  1) 

Alternatively,  if  4 TTn,R/X  » 1 the  particle  acts  as  a surface  absorber. 

4 7 

Whenever  none  of  the  above  mentioned  limiting  cases  applies,  Mie  theory 

should  be  used. 
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The  optical  constants  for  A1203  are  required  input  for  any  deter- 
mination of  the  amount  of  laser  energy  absorbed.  Unfortunately  n,  and  n, 
vary  rapidly  with  temperature,  wavelength  and  method  of  experimental 
determination.  The  values  reported  by  several  authors  * ’ are  listed 
in  Table  2„  1.  Clearly  the  experiments  based  upon  particulates  differ 
from  the  results  obtained  from  bulk  samples.  According  to  Toon  and 

g 

Pollack,  experiments  using  particulates  are  difficult  to  interpret  correctly. 

We  therefore  base  our  estimates  upon  the  bt.  Ik  sample  experiments.  The 
values  for  n2  from  Ref.  10  are  shown  in  Fig.  2.  1 for  various  temperatures 
and  wavelengths.  The  imaginary  point  of  the  index  of  refraction  is  seen 
to  increase  rapidly  with  wavelength  near  10.6^lm  for  A1  O,  at  room 
temperature.  At  higher  temperature,  however,  the  increase  with  wavelength 
is  less  significant.  In  breakdown  studies  the  particle  is  heated  from  room 
temperature  to  the  vaporization  temperature  of  4000°K.  Clearly  reliable 
values  of  the  index  of  refraction  over  the  full  temperature  range  from 
300°K  to  4000°K  are  a prerequisite  for  accurate  predictions.  In  their 
absence,the  values  of  n^  = . 8 and  n2  = . 2,  which  were  reported  in  Ref.  10 
at  T = 1773°K,  have  been  used  to  characterize  the  entire  heating  process. 

This  choice  is  logical  because  most  of  the  energy  needed  to  vaporize  the 
particulate  must  be  supplied  at  the  boiling  temperature;  therefore,  the  high 
temperature  values  of  n2  are  most  important. 

For  n2  = . 2 and  a radius  of  15  U,  the  product  of  absorption  coefficient 

and  R for  10.  6^1  radiation  is  3.  6.  Mie  theory  should  normally  be  used  in 

order  to  describe  this  regime.  Unfortunately  the  standard  Mie  theory 

code**  has  not  been  tested  for  n^  < 1,  and  results  predicted  by  it  are  not 

guaranteed  reliable.  Furthermore  the  code  provides  information  on  the  scatter- 

12 

ing  rather  than  the  local  deposition  within  the  particulate.  Instead  of 
modifying  the  code  extensively  to  overcome  these  shortcomings,  we 
resorted  to  ray  optics  to  predict  local  energy  deposition.  ^ The  Mie 
code  was  used  to  estimate  total  absorption  and  to  lend  credance  to  the  ray 
theory  prediction. 
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Fig.  2.  i Imaginary  Part  of  the  Index  of  Refraction  of  Al^O^  as  a 
Function  of  Wavelength  for  Various  Tempor&tuves. 
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Sample 


In  ray  optice  the  path  of  the  radiation  through  the  particulate  is 
determined  by  the  real  part  of  the  index  of  refraction  and  by  Snell*  s law. 


n.  sin  8 = sin  0. 
1 A 


(2.2) 


where  0 ^ is  tire  angle  between  the  ray  in  the  air  and  the  normal  to  the 
surface,  9 is  the  angle  between  the  ray  in  the  particle  and  the  normal  to 
the  surface,  and  the  index  of  refraction  of  air  has  been  set  to  one.  Typical 
rays  for  n^  1 and  nj  < 1 are  shown  in  Fig.  2.  2a  and  2.  2b  respectively. 

The  absorption  of  the  radiation  is  characterized  by  the  bulk  attenua- 

dl 

tion  coefficient  a = The  local  rate  of  deposition  is  T“*  = ol  where  I 

is  the  local  radiation  flux.  In  terms  of  the  distance  z along  the  ray  from 
the  particulate  surface  and  the  incident  laser  flux,  den  d by  I , the  local 
deposition  attributable  to  a given  ray  is  a IQ  e , Therefore  the  total 
absorption  of  the  particle  can  be  written  as 


2 R cos  9 


, -az  T 
dz  a e 1 


where  y is  the  displacement  of  the  incident  ray  from  the  ray  passing  through 
the  center  of  the  particle  as  shown  in  Fig.  2.2.  Application  of  Snell's  law 
yields 


• 2a 
sin  w 


R cos  0 = R 


_ 2 2 , 2 
R - y !nl 


The  maximum  value  of  y in  the  integration  is  R if  n^  > 1 and  R n^  if  n^  < 1. 
The  total  energy  absorbed  by  the  pail'vulate  can  then  be  expressed  as 
IoTTR  F,  where  the  efficiency  factor  F is  defined  as  1 


• ’ • V.T  1's.^'Y 
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(2.3) 


Results  obtained  by  using  Eq.  (2.  3)  with  n^  = . 8,  n^  = . 2 and 
X = 10.  6 p are  given  in  Fig.  2.  3.  The  values  predicted  by  the  Mie  code 
are  shown  for  comparison.  Three  limiting  cases  are  also  plotted;  they 
are  (1)  the  Rayleigh  limit 


> 


where  n = - i n^. 


(2)  the  surface  absorber  limit 


9 


and  (3)  the  volume  absorber  limit 


F = 


4R 

3 


a 


The  Mie  code  approaches  the  Rayleigh  limit  for  small  R,  thereby  suggesting 
that  the  code  is  reliable  even  for  < 1.  The  ray  optics  prediction  approaches 
the  surface  absorption  limit  for  large  R and  volume  absorption  limit  for 
small  R.  This  behavior  for  small  R is  the  correct  mathematical  limit  for 
ray  theory,  but  is  not  the  correct  limit  for  the  physical  phenomenon  since 
the  requirement  X « R is  no  longer  met.  In  the  regime  of  interest,  R = 15 
ray  theory  predicts  61%  absorption  whereas  Mie  theory  predicts  87%  absorption 
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ch  is  Absorbed  by  a S 


The  value  from  Mie  theory  will  be  used  for  calculations  which  require  only 
the  total  energy  absorbed.  For  calculations  concerned  with  local  energy 
deposition  the  ray  theory  predictions  will  be  used  to  suggest  the  maximum 
energy  deposition  rates  within  the  particle.  The  assumption  we  make  in 
this  approach  is  that  the  enhanced  absorption  of  Mie  theory  comes  from 
stronger  absorption  in  regions  where  ray  theory  does  r><-t  predict  much 
absorption,  whereas  the  maximum  absorption  rates  whxcn  generally  prevail 
are  unaltered. 

The  discrepancy  between  the  two  theories  suggests  that  calculations 
accurate  to  better  than  30%  require  a complete  analysis  using  full  Mie 
theory  to  study  deposition.  Such  an  ambitious  project  should  not  be  under- 
taken,  however,  without  also  understanding  how  to  incorporate  the  effects 
of  both  spatial  variations  of  the  index  of  refraction  (caused  by  non-uniform 
heating)  and  non-uniform  surface  recession  (caused  by  vaporization).  (Of 
course  many  other  factors,  such  as  particle  motion,  which  are  neglected 
in  this  report,  are  also  important  for  more  accurate  calculations. ) 

Prediction  of  the  time  when  vaporization  commences  requires 

information  about  local  energy  deposition.  We  can  qualitatively  study 

local  deposition  via  ray  theory.  In  Fig.  2.  4 several  rays  are  shown 

striking  an  AlgO^  particle.  The  path  of  the  ray  directed  towards  the  center 

is  not  changed,  therefore  the  energy  deposition  along  its  path  is  proportional 

to  ae_CtZ  where  z is  the  distance  along  the  ray  from  the  front 

surface.  Two  competing  processes  influence  the  deposition  rate  at  other 

interior  point 8.  The  convex  shape  of  the  interface  causes  rays  to  diverge 

after  passing  into  the  particulate,  hence  the  intensity  in  the  interior  drops 
• ctz 

faster  than  e due  to  the  geometric  spreading.  On  the  other  hand,  the 
initial  effect  of  the  change  in  direction  on  passing  through  the  interface  is 
to  squeeze  the  rays  together,  thereby  creating  a more  intense  beam.  There- 
fore, at  the  surface,  the  intensity,  and  hence  the  absorption  per  unit  volume,  is 
enhanced.  The  numbers  in  Fig.  2.4  associated  with  each  ray  in  the  relative 
intensity  immediately  inside  the  particulate.  The  numbers  along  the  central 
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ray  represent  the  intensity  at  interior  points  if  defocussing  is  neglected. 

It  is  evident  that  "hot  spots"  exist  where  energy  is  absorbed  much  more 

quickly  than  over  the  rest  of  the  surface.  The  bulk  of  the  particle  surface, 

however,  absorbs  at  a rate  only  minimally  above  the  absorption  rate  for  the 

j central  ray.  Therefore  bulk  heating  is  characterised  by  a Iq. 

The  final  ingredients  needed  to  predict  vaporisation  time  scales 

are  the  thermodynamic  properties  of  A1203  in  its  various  forms.  The 

JANAF15  tables  were  used  as  sources  of  melting  temperature,  specific 

heat,  free  energy  and  enthalpy  of  formation  for  solid  A1203,  liquid  A1203 

and  the  vapor  components  AlO,  A102,  A120,  A1202,  02>  O^,  AlO  , A^O  , 

A1  O O",  O ~ and  Al.  O ".  The  energy  levels  and  degeneracies  of 

2 2.,  2 2 2 + ++  + 

Moore  were  used  for  the  atomic  species  - Al,  Al  , Al  , O and  O • 

The  vapor  pressure  above  the  Al^  liquid  surface  at  any  temperature  T 

was  determined  by  an  equilibrium  code  calculation  which  solves  the 

equations  relating  species  concentrations  in  thermodynamic  equilibrium. 

A plot  of  vapor  pressure  versus  T is  given  in  Fig.  2.  5,  along  with  the 

enthalpy  of  the  vapor.  The  boiling  point  is  thus  determined  to  be  4000  K 

and  the  heat  of  vaporization  at  4000°K  is  found  to  be  2.  02  x 10  j/mole. 

It  can  be  seen  from  Fig.  2.  5 that  the  enthalpy  per  unit  mass  does  not 

change  appreciably  over  the  temperature  range  displayed.  A summary 

of  useful  thermodynamic  properties  of  Al-O  is  given  in  Table  2.  2. 

6-3  17 

Thermal  conductivity  and  density  were  taken  from  the  AIP  handbook. 

The  time  T at  which  the  surface  begins  to  vaporize  is  given  by 
v 

t = p Ah/a  I (2. 4) 

v o 
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TABLE  2.  2 


Proparties  of  A^O^  Used  in  Calculation 


Quantity 

Symbol 

Value 

Ref. 

Density 

P 

3.  8 g/cc 

17 

Specific  heat  (average) 

c 

p 

1.4  .T/gK 

15 

Atomic  weight 

A 

101.96 

15 

Thermal  conductivity 

(average) 

X 

c 

. 1 W / cm  K 

17 

Melting  temperature 

T 

m 

2315  K 

15 

Normal  boiling  point 

T 

V 

400  K 

equilibrium  code 

Heat  of  fusion 

h 

m 

1.  19  x 103  J/mole 
1.67xl03J/g 
4.  43  x 103  J/cm3 

* 

15 

Heat  of  vaporisation 

h 

V 

2.  02  x 10^  J/mole 
1.  98  x 104  J/g 
7.  53xl04  J/cm3 

equilibrium  code 

Enthalpy  difference  between 
gas  at  4000°K  and  1 atm 
and  solid  at  300°K 

Ah 

! 

2.  63  x 10^  J/mole  I 
2.  58  x 104  J/g  | 
9.8xl04J/cm3  ] 

equilibrium  code 
and 
15 

Real  part  cf  index  of  refrac- 
tion (average) 

nl 

.8 

10 

Imaginary  part  of  index  of 
refraction 

n2 

. 2 

10 

Bulk  absorption  coeffient, 
X = 10.  6u 

a 

2370  cm  * 
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where  p ie  the  solid  density  end  Ah  is  the  energy  required  to  change 
Al^O^  from  a solid  at  room  temperature  to  a vapor  at  the  vapor  tempera- 
ture. For  the  purposes  of  this  calculation  Al^Og  is  assumed  to  vaporise 
directly  into  the  equilibrium  vapor  constituents.  It  is  shown  in  Fig.  2.  5 
that  the  enthalpy  of  the  vapor  does  not  change  appreciably  as  the  vapor 
temperature  is  raised.  Therefore  a constant  value  of  Ah  = 2.63x10^ 

4 

J/mole  = 2.  58  x 10  J/g  is  used  to  calculate  T • This  corresponds  to  a 
vapor  pressure  of  one  atmosphere  and  a vapor  temperature  of  4000°K. 

We  find  for  X = 10.  6|J  that  the  fluence  required  to  vaporize  the  surface  is 

!o  Tv  » 41*4  J/cm2  , (2.5) 


where  I is  the  incident  laser  flux.  In  the  range  of  intensity  of  interest, 
bulk  vaporization  commences  from  5 to  14  microseconds  after  start  of 

2 

the  pulse.  As  an  aside  it  should  be  noted  that  the  requirement  of  41  J/cm 

to  cause  bulk  vaporization  does  not  conflict  with  the  observed  breakdowns 

2 

shown  in  Fig.  1.  1 which  occur  with  only  10  J/cm  fluence.  These  observa- 
tions are  made  at  high  intensities  where  bulk  vaporization  is  not  required 

to  initiate  breakdown.  Indeed,  the  "hot  spots"  mentioned  earlier  vaporize 
2 

before  10  J/cm  is  received  and  their  vaporization  may  be  sufficient  to 
induce  breakdown  at  high  intensity. 

The  calculation  neglected  thermal  conduction  within  the  particulate 
and  conduction  losses  at  the  surface  to  the  air.  The  losses  to  the  air  are 
minimal.  As  will  be  shown  in  the  Section  3.  5 the  flux  at  the  surface  into 
the  air  at  time  t is  given  by 


R 

1 + >J  fT  Kt 


9 


(2.6) 


^ *2 

where  X is  the  thermal  conductivity  of  air  (5.  85  x 10  W/cmK  at  4000°K 
and  one  atmosphere),  AT  is  the  temperature  difference  between  the  surface 
and  the  initial  air  temperature  and  K is  the  thermal  diffusivity  of  air 
2 0 

(32.8  cm  /sec  at  4000  K and  one  atmosphere).  The  maximum  flux  occurs 

when  the  liquid  reaches  the  vaporization  temperature,  but  before  any  of 

o 5 2 

the  heat  of  vaporization  is  added.  For  T = 4000  K,  I = 5 x 10  W/cm  , 

° 4 

AT  = 3700  K,  this  time  is  given  by  t ~ 9.  6/ 1 = 2 Msec.  Then  q = 1.  6 x 10 

2 

W/cm  is  the  maximum  flux  expected  to  be  lost  through  conduction  to  the 

air.  However  the  rate  of  energy  deposition  immediately  below  the  surface 
10  3 

is  al  = 1.  2 x 10  W/cm".  Clearly  the  conduction  losses  can  be  supplied 

by  the  energy  deposited  in  a very  thin  surface  layer. 

A potentially  more  serious  effect  is  thermal  conduction  within  the 

particle.  The  thermal  diffueivity  for  A^O^  is  found  from  Table  2.  2 to  be 
-2  2 

1.8  x 10  cm  /sec.  If  thermal  conduction  is  ineffectual,  the  temperature 
profile  at  time  t is  given  by 


I, 

t a e 

T U)  = -2— 

PCP 

where  z is  the  distance  from  the  surface.  For  this  profile,  the  local  source 
term  due  to  thermal  conduction  is 


■Six  = a3  I t e-°z  K - 
d*2 


TC 


which  is  to  be  compared  to  the  source  term  due  to  laser  absorption 

-GLX 

S,  = a I e * Their  ratio  is 
L o 


STC/SL=ttKt 
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As  before  the  time  over  which  the  sources  are  to  be  compared  is  the  time 

6 2 

required  to  reach  the  vaporization  temperature.  For  I = 3 x 10  W/crn 

-6  0 
the  time  is  3 x 10  sec  which  leads  to  a ratio  of  . 3.  Therefore,  even 

at  its  maximum  value,  thermal  conduction  is  small  compared  to  laser 

heating. 

Having  established  the  shortest  time  for  bulk  vaporization  to  occur, 
we  now  characterize  the  time  to  vaporize  a given  fraction  (3  of  the  particle. 
In  developing  a simple  model  for  this  process  we  assume  that  all  the 
absorbed  energy  is  effective  in  vaporizing  the  particulate.  Under  the 
further  assumption  that  the  absorption  cross  section  does  not  change 
significantly  during  vaporization  we  find  the  relationship 


teI  ttR2F  = P p A h ■■  "n,  "■ 

P O 3 


(2.7) 


where,  as  before  F is  the  efficiency  factor  and  Ah  is  the  energy  required. 
For  a 1 .5  p.  particle  of  Al  O with  F = . 87,  the  time  to  change  a fraction  (3 
of  the  particle  to  vapor  at  one  atmosphere  and  4000  K is  given  by 


I t-  = 225  j3  J/cm 
o p 


(2.8) 


Of  course  for  small  {3  the  prediction  is  obviously  Incorrect  - vaporization  has 
not  begun.  In  the  limit  of  small  |9  the  correct  time  is  the  sum  of  tg  and  t^. 
The  approximation  should  be  reliable  as  |3  approaches  one,  however,  since 
then  all  the  absorbed  energy  then  contributes  to  the  vaporization.  The  time 
is  predicted  by  Eq.  (2.  7)  is  tabulated  in  Table  2.  3 for  three  values  of  Iq 
and  three  values  of  j3.  Clearly  a major  portion  of  the  observed  breakdown 
time  is  spent  producing  the  vapor.  It  should  be  kept  in  mind  that  the  times 
predicted  by  Eq,  (2.  7)  are  only  rough  estimates.  It  is  not  known  how  the 
energy  absorbed  changes  once  vaporization  begins  to  change  the  geometry 
of  the  particulate. 
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The  second  stage  of  the  breakdown  process  presupposes  a pure  cloud 

of  the  vapor  exists.  However,  a pure  cloud  cannot  be  formed  unless  the 

vaporization  rate  is  much  faster  than  the  rate  at  which  the  vapor  can  diffuse 

into  the  air.  It  is  common  to  assume  that  a pure  cloud  will  be  formed  if 

the  pressure  of  the  vapor  at  the  particle  surface  is  above  two  or  three 

atmospheres.  An  estimate  of  the  pressure  of  the  vapor  at  the  interface 

can  be  made  as  follows:  The  mass  flux  is  determined  by  the  heating  rate. 

In  the  limit  of  one-dimensional  steady  state  recession  the  mass  flux  per 

unit  area  p u is  related  to  the  laser  intensity,  which  is  assumed  to  be 
s s 

entirely  absorbed,  by 


p u Ah  = I , 
s s o 

where,  as  before.  Ah  is  the  enthalpy  required  to  change  a solid  to  a vapor 

at  the  required  conditions.  The  temperature  and  vapor  pressure  are  related 

as  shown  in  Fig.  2,  5,  For  the  present  estimates,  a constant  value  of 

£ 

Ah  = 2.63  x 10  J/mole  can  be  used.  Under  the  assumption  that  there  is 
massive  vaporization,  the  velocity  of  the  vapor,  denoted  by  u^,  assumes 
its  limiting  value,  the  sonic  velocity,  which  defined  by 

uv  = /yrt7  , 

where  R is  the  ideal  gas  constant  and  Y is  the  ratio  c^/c^.  Furthermore, 
using  both  mass  continuity  at  the  interface, 

p u = p u , 

8 S rV  V 

where  p^  is  the  vapor  density,  and  the  ideal  gas  law, 

P = Pv  RT  * 

obtain  the  approximate  relationship 
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or 


(2,9) 


p [atm]  = 8.  2 x 10  ^ I [W/cm^]  . 


(2.  10) 


Equation  (2.9;  5 s expected  to  predict  reasonable  values  of  the  pres- 
sure whenever  steady  state  recession  is  approached.  Immediately  after 
vaporization  commences,  the  pressure  is  higher  than  the  steady  ate  value, 
but  for  the  time  scale  required  to  obtain  appreciable  vaporization  the  steady 
state  pressure  is  representative  of  the  average  pressure. 

Several  values  of  p as  determined  by  Eq.  (2.  10)  are  given  in 
Table  2.4.  The  equation  breaks  down  for  low  pressures  because  u^  is 
no  longer  the  sonic  velocity.  It  should  hold,  however,  for  the  pressures 
of  interest,  that  is,  for  pressures  greater  than  two  or  three  atmospheres. 

From  Eq.  (2. 10)  then,  a reasonable  threshold  intensity  for  production  of 

6 2 

a pure  vapor  is  3 x 10  W/cm  . This  criteria  agrees  well  with  the  observed 
threshold,  which  suggests  that  a pure  vapor  cloud  is  indeed  needed. 

The  pressure  calculation  does  not  involve  the  particle  radius;  it  only 
demands  that  the  particle  be  large  enough  it  o justify  the  assumption  of 
steady  state  recession.  Therefore  a lower  limit  has  been  calculated  which 
is  valid  for  the  production  of  a pure  A1  O vapor  over  even  a semi  infinite 

C*  J 

plane.  If  breakdown  does  indeed  require  a pure  cloud,  then  large  particles 
should  have  the  same  intensity  threshold  as  the  30  ^ diameter  particles  of 

a:2o3. 

To  summarize,  the  rough  estimates  of  the  time  to  the  start  of  vapor- 
ization, the  time  to  vaporize  an  arbitrary  fraction  of  the  particle,  and  the 
vapor  pressure  above  the  interface  are  consistent  with  experimental  observa- 
tions of  the  delay  time.  Furthermore,  it  has  been  shown  that  the  threshold  for 
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TAB  LE  2.  4 


Pressure  as  a Function  of  10.  6 p.  Laser  Intensity  for 
- Assuming  Steady  State  Recession  and  Rapid  Vaporization 


[W /cm^] 

p [atm] 

T 

vapor 

io6 

. 82 

3960 

2 x 106 

1.6 

4110 

3 x 106 

2.  5 

4210 

4 x 106 

3.  3 

4280 

5 x 106 

4.  1 

4340 

6 x 106 

4.9 

4390 

7 x 106 

5.  7 

4430 

8 x 106 

6.6 

4470 
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production  of  pure  Al^O^  vapor  is  about  3 x 10  W/cm  for  10.  6 M radiation. 
Finally  the  time  required  to  produce  the  vapor  cloud  is  a significant  portion 
of  the  delay  time  observed  before  breakdown. 

The  final  subject  to  be  discussed  in  this  section  is  how  these  predic- 
tions for  vapor  production  change  as  the  laser  wavelength  and  particle 
composition  are  altered.  For  the  calculation  of  t the  important  quantities 
are  the  energy  required  to  vaporize  a unit  volume,  p Ah,  and  the  energy 
deposition  rate  per  unit  volume,  denoted  by  a I.  It  is  important  that 
the  maximum  value  of  al  be  used.  For  1 the  maximum  occurs  at  the 

front  surface,  but  for  ot  R « 1 the  maximum  value  may  occur  at 
the  back  of  the  particle  because  of  focusing  effects.  SiO  particles  may  show 

Ct 

this  effect.  The  time  to  vaporize  a given  fraction  0 of  the  particle  is  once 
again  given  by  Eq.  (2.7).  Only  the  composition  of  the  particulate  affects 
p Ah,  but  both  the  composition  and  the  wavelength  affect  F.  The  pressure 
of  the  vapor  is  related  only  to  the  composition,  under  the  assumption 
that  the  recession  rate  is  steady.  This  assumption,  however,  is  valid 
only  if  the  laser  energy  is  absorbed  reasonably  close  to  the  surface; 
therefore  the  applicability  of  Eq.  (2.9)  depends  upon  laser  wavelength 
as  well  as  composition.  As  a final  comment,  it  should  be  noted  that 
Al  O has  a small  imaginary  index  of  refraction  at  short  wavelengths. 

In  fact,  at  3.  8pi,the  bulk  absorption  coefficient  is  of  the  order  of  . 3 cm 
rather  than  2370  cm  \ The  particle  is  effectively  transparent  and  does 
not  vaporize.  Long  pulse  Al  O induced  breakdown  at  3.  8 |i  will  not  occur 
at  the  intensities  studied  here,  although  it  may  still  occur  at  these  intensities 
for  other  particulates. 


! 


3.  VAPOR  CLOUD  HEATING 


The  second  stage  of  the  long  pulse  breakdown  process  is  the  heating 
of  the  vapor  cloud  by  the  laser.  In  this  section  two  simple  models  are  used 
to  estimate  both  the  time  to  heat  the  vapor  to  20,  000°K  and  the  threshold 
intensity  below  which  the  vapor  will  not  heat.  One  model  is  appropriate 
for  short  times;  it  assumes  that  the  vapor  in  the  center  of  the  cloud  is 
effectively  insulated  from  the  cold  air  by  the  outer  layers  of  vapor.  The 
other  model  is  appropriate  for  longer  times  and  slow  heating.  Thermal 
conduction  is  able  to  establish  a temperature  profile  wherein  the  losses 
are  effectively  distributed  over  the  complete  vapor  cloud.  In  both  models 
the  heating  is  assumed  to  occur  slowly  enough  that  the  pressure  remains 
constant  and  the  composition  of  the  vapor  cloud  is  determined  by  thermo- 
dynamic equilibrium.  The  constant  pressure  assumption  is  reasonable 
since  the  heating  is  known  to  take  several  microseconds,  whereas  the 
characteristic  time  for  equilibriating  pressure  is  of  the  order  of  . 3 micro- 
seconds. The  chemical  equilibrium  assumption,  at  least  in  the  beginning 
stages,  involves  the  rate  of  dissociation  of  molecules.  Characteristic 
times  for  establishing  chemical  equilibrium  at  4,  000°K  are  also  expected 
to  be  less  than  a microsecond. 

3.  1 Initial  Conditions 

Before  the  models  can  be  exercised  the  initial  conditions  must  be 
determined.  Since  the  pressure  equilibrates  with  the  surrounding  air  within 
. 3 microseconds,  the  vapor  cloud  starts  (and  remains)  at  one  atmosphere. 
The  initial  temperature  is  then  expected  to  be  near  4000°K  - otherwise  the 
vapor  should  liquify.  A more  accurate  estimate  of  the  initial  temperature 
can  be  calculated  as  follows:  Given  a value  of  the  incident  laser  intensity  I , 


-29- 


3.  VAPOR  CLOUD  HEATING 


The  second  stage  of  the  long  pulse  breakdown  process  is  the  heating 
of  the  vapor  cloud  by  the  laser.  In  this  section  two  simple  models  are  used 
to  estimate  both  the  time  to  heat  the  vapor  to  20,  000°K  and  the  threshold 
intensity  below  which  the  vapor  will  not  heat.  One  model  is  appropriate 
for  short  times;  it  assumes  that  the  vapor  in  the  center  of  the  cloud  is 
effectively  insulated  from  the  cold  air  by  the  outer  layers  of  vapor.  The 
other  model  is  appropriate  for  longer  times  and  slow  heating.  Thermal 
conduction  is  able  to  establish  a temperature  profile  wherein  the  losses 
are  effectively  distributed  over  the  complete  vapor  cloud.  In  both  models 
the  heating  is  assumed  to  occur  slowly  enough  that  the  pressure  remains 
constant  and  the  composition  of  the  vapor  cloud  is  determined  by  thermo- 
dynamic equilibrium.  The  constant  pressure  assumption  is  reasonable 
since  the  heating  is  known  to  take  several  microseconds,  whereas  the 
characteristic  time  for  equilibriating  pressure  is  of  the  order  of  . 3 micro- 
seconds. The  chemical  equilibrium  assumption,  at  least  in  the  beginning 
stages,  involves  the  rate  of  dissociation  of  molecules.  Characteristic 
times  for  establishing  chemical  equilibrium  at  4,  000°K  are  also  expected 
to  be  less  than  a microsecond. 

3.  I Initial  Conditions 

Before  the  models  can  be  exercised  the  initial  conditions  must  be 
determined.  Since  the  pressure  equilibrates  with  the  surrounding  air  within 
. 3 microseconds,  the  vapor  cloud  starts  (and  remains)  at  one  atmosphere. 
The  initial  temperature  is  then  expected  to  be  near  4000°K  - otherwise  the 
vapor  should  liquify.  A more  accurate  estimate  of  the  initial  temperature 
can  be  calculated  as  follows:  Given  a value  of  the  incident  laser  intensity  I , 
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the  vapor  temperature  and  the  pressure  of  the  vapor  at  the  surface 

of  the  particulate  are  found  from  Fig.  2.  5 and  Eq.  (2.  10)  respectively. 

The  vapor  expands  to  one  atmosphere  in  a fraction  of  a microsecond. 

Therefore,  the  laser  heating  during  the  expansion  should  be  negligible 

and  the  initial  temperature  Tg  of  the  vapor  cloud  corresponds  to  isentropic 

expansion  of  the  vapor  from  p^and  T^to  one  atmosphere  andT^.  A plot 

of  T^  versus  p^  for  Al^O^,  as  determined  by  the  equilibrium  code,  is 

given  in  Fig.  3.  1.  Of  particular  interest  are  the  values  3940°K  and  3870°K 

which  are  reached  when  the  laser  intensities  are  3 x 10  W/crn  and 
6 2 

8 x 10  W/cm  respectively.  These  values  straddle  the  re  sting  range 

of  intensity.  Since  the  values  of  Tg  lie  so  close  to  the  normal  boiling 
point  the  calculations  performed  in  this  section  are  made  with  a constant 
initial  temperature  of  4000°K.  For  completeness,  however,  the  various 
physical  parameters  are  evaluated  for  temperatures  as  low  as  3500°K. 

The  last  parameter  to  be  determined  is  the  initial  vapor  cloud  size.  A 
fraction  0 of  a Al^O  particle  which  is  15  p in  radius  produces  a vapor 

o * 1/3 

cloud  at  4000  K and  one  atmosphere  which  has  a radius  of  3 S 2(3 ) microns. 
3.  2 Physical  Parameters  Required  by  the  Models 

The  local  equation  governing  the  temperature  history  of  the  vapor 
is 

p (T)  cp  (T)~=  <T)Iq  -V*  S-  V«  \c(T)7T  (3.  1) 

« 

where  p is  the  density,  cp  the  specific  heat,  the  effective  laser  absorp- 
tion coefficient,  S the  radiative  flux  produced  by  the  vapor  and  X the 
coefficient  of  thermal  conductivity.  Approximations  for  V * S and  7*  V T, 
which  are  made  later  in  this  section,  depend  upon  the  vapor  cloud  radius 
R (T).  The  use  of  the  incident  laser  intensity  Iq  throughout  the  entire 
volume  assumes  that  the  cloud  is  transparent  to  laser  radiation.  The 
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Fig.  3.  1 The  Relationship  Between  the  Pressv'e  of  a AI2O3  Vapor  Cloud  in 
Equilibrium  with  Liquid  AI2O3  and  the  Temperature  of  the  Vapor 
after  Isentropic  Expansion  to  One  Atmosphere. 
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values  of  K ^ and  R determined  in  the  following  subsections  justify  this 
simplification.  In  the  subsequent  subsections,  the  physical  parameters  p, 

i _ _ 

c , R(T),  K (T)  and  X (T)  are  discussed,  and  the  models  for  7 * S and 
p l c 

T • X ?T  are  described, 
c 

3.  3 Parameters  Determined  bv  the  Equilibrium  Code 

The  equilibrium  code  gives  the  density,  enthalpy  and  specific  heat 
of  the  vapor  in  addition  to  the  concentrations  of  the  various  species.  In 
Fig.  3.  2 we  plot  the  proc  uct  p c^.  The  large  values  at  low  temperature 
indicate  that  a large  fraction  of  the  total  energy  must  be  supplied  in  the 
initial  heating.  In  Figs.  3.  3 and  3.4  the  mole  fraction  of  the  various 
species  are  shown.  Figure  3.  3 gives  the  composition  for  low  temperature 
where  many  molecular  species  are  present  and  rapid  variations  in  the 
composition  occur  within  a few  hundred  degrees.  Figure  3.  4 displays  the 
relative  concentration  of  ions,  electrons  and  atoms  over  the  complete 
range  of  temperature  of  interest.  The  density  of  the  vapor  can  be  used  to 
determine  the  radius  of  the  vapor  cloud.  In  Fig.  3.  5 the  radius  is  given 
as  a function  of  temperature  for  three  values  of  3. 

3.4  Laser  Absorption  Coefficient 

Two  types  of  transitions  can  absorb  laser  radiation:  free-free  transi- 
tions of  electrons  during  collisions  with  heavy  particles,  and  molecular 
transitions  between  states  of  different  vibrational  and  rotational  quantum 
numbers.  The  electron  transitions  - called  inverse  Bremsstrahlung  - absorb 
for  any  wavelength  of  laser  and  their  absorption  cross  section  scales  roughly 
as  wavelength  squared.  Molecular  transitions,  on  the  other  hand,  occur 
over  limited  frequency  bands,  and  a change  of  laser  wavelength  may  strongly 
enhance  the  absorption  or  eliminate  it  entirely. 

The  inverse  Bremsstrahlung  absorption  coefficient  is  given  by 


* ' ..  -1.43883  u)/T. 
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where  uu  is  the  wave  number  of  the  laser  radiation  in  cm  \ T is  the 

temperature  in  °K,  n is  the  concentration  of  electrons,  n.  the  concentre- 

th  ® ^ 

tion  of  the  i species  and  Q.  the  average  cross  section  for  absorption  of 

a photon  of  wavenumber  U)  by  an  electron  colliding  with  a particle  of  species  i. 

For  charged  species  the  cross  section  can  be  approximated  by  the  product 

IS 

of  the  Gaunt  factor  g,  given  by  Karzas  and  Latter,  and  the  semi  classical 

19 

formula  of  Kramer: 


0 4 — / 2tt 

Qi  " 3 *•  ^3  mg  kT 


1/2 


6 7 2 
e Z, 
i 


4 . 3 

c m.  h w 
e p 


(3.3) 


where  uu  is  the  wave  number  of  the  radiation,  m is  the  mass  of  the  electron, 

e 

hp  is  Planck's  constant,  k is  Boltzmann's  constant,  e is  the  charge  of  the 
electron,  c is  the  velocity  of  light  and  is  the  charge  of  the  particle  in 
units  of  electron  charge. 

The  cross  section  of  the  neutral  species  is  not  accurately  known 
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although  typical  values  are  of  the  order  of  10  cm  . We  used  the  following 

estimates  of  Q for  the  neutral  species.  For  Al,  a highly  polarizable  atom, 

20 

the  expression  of  Hyman  and  Kivei  and  the  polarizability  given  by  Teachout 
21 

and  Pack  were  used  to  nredict  the  cross  section 


QA1  = .64  (T/103)  10"36 


(3.4) 


The  crons  sections  for  O and  O2  were  deduced  from  the  experimental  values 

22 

given  by  Taylor  and  Caledonia.  All  others  were  guessed  by  using  the 

predictions  for  Al,  O and  O as  guides.  The  cross  sections  used  herein  are 

listed  in  Table  3.  1.  The  prediction  of  neutral  and  charge  inverse  Brem- 

sstrahlung  absorption  coefficients  for  I 0.  6 p.  radiation  and  equilibrium  Al  0„ 

vapor  at  one  atmosphere  are  plotted  in  Fig.  3.  6.  The  total  inverse  Brem- 

sstrahiung  absorption  coefficient  is  shown  as  well.  The  rapid  increase  of 

K*  with  the  temperature  suggests  that  the  initial  heating  of  the  vapor  accounts 
-Ll 

for  most  of  the  time  consumed  in  the  vapor  heating  stage.  Electron-neutral 
inverse  Bremsstrah lung  dominates  in  the  low  temperature  regime. 
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TABLE  3.  1 


Species 

Estimated  Q for  10.  6 n Radiation 

[cm5] 

A1 

. 64  (T/103)  10'36 

O 

3 x 10" 36 

°2 

1.2xl0‘35 

AlO 

6 x 10-36 

A12C 

1.2  x 10"35 

aio2 

1.  2 x 10"35 

A12°2 

-35 

1.  2 x 10 
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The  other  mechanism  whereby  10.  6|i  radiation  can  be  absorbed  in 
molecular  absorption.  Among  the  molecular  species  in  Al^O^  vapor  only 
AlO  has  its  vibration- rotation  band  near  10.  6|i.  In  order  to  estimate  the 
absorption  accurately  the  precise  wave  number  of  the  laser  line  is  required. 

The  laser  used  in  breakdown  experiments  lases  mainly  on  the  P(20)  line  of 

23 

CO  , although  it  often  fires  on  the  P(18)  and  P(22)  line  also.  The  wave 

2 24 

numbers  of  these  lines,  as  given  by  McClatchy  and  Selby  are  listed  in 

Table  3.  2.  For  the  temperature  and  pressure  of  interest  the  lines  of  the 

AlO  vibration- rotation  transitions  are  narrow;  therefore  the  absorption 

coefficient  must  be  calculated  by  summing  the  contributions  of  individual 

lines  rather  than  by  using  an  overlapping  line  model.  We  used  the  following 

2 + 

method  to  compute  the  absorption.  The  energy  levels  of  the  AlO  X £ state 
are  given  by 


e =v(u)  - u)x)-v2(dx  + B J (J+  1)  - a (v  +~)  J (J  + 1) 
v,  J eee  eee  e2 


(3.5) 
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C u 6 


for  v i O and  J i 0 . 


The  zero  of  energy  is  chosen  to  be  € , and  the  values  of  the  constants 
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as  reported  by  Sulzman  are  listed  in  Table  3.  3.  The  R branch  transition 

g 

(Av  = 1,  AJ  = 1)  from  the  v,  J state  has  wave  number  ujv  j which  is  given  by 


J*  = uj  - 2U)  x (v  + 1)  + 2 (V  + 1)  [B 
, J e e e < 

-4TD  +(v+|)p  ][j+i]3-e  (J+1)2 
e 4e  e 


-(v+-;)a  ] - a [J+l][J+2] 
2 e e 

(J  + 2)2  for  V 2 0,  JJs  0 


(3.6) 
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whereas  the  P branch  transition  (Av  = 1,  AJ  = - 1)  has  wave  number  oj  1 

v,  J 

given  by 


Ui.  f = (JU  - 2uj  x (v  + 1)-2J  [B  - (V  +r)a  ] - J(J-  l)a  +4j3[d  +(v+4)P  1 

V * J 6 G G 0 £0  0 0 2 e 

(3.7) 

for  v * 0 and  J a 1 . 


-J2  (J  + l)2  0 


The  Q branch  transition  (Av  = 1*  AT  = 0 ) is  not  present  fora 
E state.  In  terms  of  the  integrated  cross  section  for  a given  R branch 
transition,  which  we  denote  by  f^  ^ , the  integrated  absorpti< 

coefficient  for  the  transition  is 


ion 
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line 
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-1.43883  s t/T 
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where  Z is  the  partition  function  for  AlO,  including  electronic, 

01^  V > J 

vibrational,  and  rotational  states,  2 is  the  electronic  degeneracy  factor 

2 + 

of  the  ground  state  X Z , n^Q  is  the  concentration  of  aluminum  oxide  mole- 
cule, e is  the  charge  of  the  electron,  m is  the  mass  of  the  electron  and 
c is  the  velocity  of  light.  An  analogous  formula  holds  for  the  P branch 
transitions. 

The  detailed  line  profile  L((u)  must  be  used  to  determine  the 
absorption  coefficient  at  the  laser  wavenumber.  Two  mechanisms  broaden 
the  molecular  rotation  lines:  (1)  Doppler  broadening  which  has  a half-half 
width  given  by 
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where  u)  ie  the  wave  number  of  the  line,  M is  the  mass  of  AlO  and  A is 
o 

its  atomic  weight;  and  (2)  collision  broadening  which  has  a half-half  width 

b . The  parameter  b , for  most  molecules,  is  of  the  order  of  . 08  cm’^ 
c c 

at  300  K and  1 atmosphere,  and  it  scales  roughly  as  1/,/F  at  constant 

26 

pressure.  In  the  regime  considered  here,  collision  broadening  is 

5 times  larger  than  Doppler  broadening,  so  that  the  Voigt  line  profile 

27 

closely  ressembles  a Lorentz  profile  with  an  effective  half-half  width 
given  by 

b = b_  + b (3.10) 

c D c 

The  profile  function  is  then  given  by 


L(UJ)=._ — ^rr (3.U) 

b4  + On  - U)  ) 

O 

The  total  effective  absorption  coefficient  for  a laser  line  is  determined 
by  summing  over  the  contributions  from  all  lines  and  by  then  multiplying 
by  (1  - e or(jer  inciude  simulated  emission.  Applying 

this  procedure  to  Eq.  (3,  11)  and  (3.  8)  we  find 
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R P 

where  b _ and  b T are  the  half-half  widths  of  the  R and  P transitions 
v* j J 

respectively,  from  the  v,  J state. 
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The  final  ingredients  needed  are  the  f numbers  of  all  the  transitions. 
Fortunately  the  f numbers  of  different  vibrational- rotational  transitions 
are  related  by  simple  formulas.  The  quantum  mechanical  formula  for  the 
f number  is 
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where  M is  the  electric  dipole  moment  of  the  transition,  m is  the  mass 
e e 

of  the  electron,  c is  the  velocity  of  light,  h is  Planck's  constant,  e is  the 

P 

charge  of  the  electron  and  uj  is  the  wave  number  of  the  transition.  Only 
u)  and  change  for  different  transitions.  For  a molecular  system  with 
the  energy  levels  approximated  by  an  harmonic  oscillator,  the  only 
transition  moment  is  between  adjacent  vibrational  states,  Av  = 1. 
Furthermore  the  strength  of  the  transition  moment  is  proportional  to 
(v  + 1).  (Hereafter  we  call  this  the  harmonic  scaling  law).  Finally  the 


relative  strength  of  the  two  allowed  angular  momentum  transitions  are  given 

J + 1 _ J „ 

by  the  Honl- London  factors:  ^ for  R transitions,  +' ^ for  P transi- 


tions. 


Therefore,  if  the  f number  for  one  transition,  say  theR  transition 
for  v = 0,  J = 0 is  known,  the  f numbers  of  all  other  transitions  are  given  by 
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The  value  of  f has  been  theoretically  predicted  by  Michels.  One 

model  he  used  predicted  f = 8.  58  x 10  ^ for  v = 0,  and  transitions  from 

high  v states  were  found  to  be  close  to  (v  + 1)  f.  Furthermore,  Av  * 1 

transitions  were  surpressed.  Unfortunately  another  model  led  to  the  value 
-5 

f = 1.25  x 10  . However,  for  the  most  important  transitions  in  our 

calculations,  the  transitions  from  the  = 5 state,  strong  anharmonic  effects 

-5 

in  the  second  method  of  calculation  reduce  the  f number  to  5.  15  x 10 

which  agrees  well  the  value  6 x 8.  58  x 10  ^ which  is  given  by  the  first 

model  and  the  harmonic  scaling  law.  Therefore,  the  f number  of 

8.  58  x 10  ^ and  the  harmonic  scaling  law  were  used  in  our  calculations. 

29 

We  assumed  that  the  value  quoted  is  correct  for  the  f transition.  and 

0,  0 

that  Eqs.  (3.  14)  define  all  the  other  required  f numbers. 

Calculations  of  the  absorption  by  vapor  for  the  P(18),  P(20) 

and  P(22)  lines  of  CO  were  performed  by  using  Eqs.  (3.  12)  and  (3.  14), 

-6  ^ 

f_  = 8.  58  x 10  and  b = . 08*/300/T  . The  results  are  shown  in  Fig.  3.7. 
0,0  c 

The  dominant  laser  line,  P(20),  is  seen  to  be  also  the  most  strongly  absorbed; 
hence  breakdown  most  likely  occurs  when  the  laser  lases  on  the  P(20)  line. 

All  subsequent  calculations  therefore  assume  the  laser  fires  on  the  P(20) 
line.  The  total  absorption  coefficient  is  plotted  in  Fig.  3.  8.  It  can  be 
seen  that  molecular  absorption  in  AlO  exceeds  inverse  Bremsstrahlung 
absorption  in  the  crucial  low  temperature  regime. 

The  molecular  absorption  coefficient  is  very  sensitive  to  two  factors  - 
(1)  the  position  of  the  laser  line  and  (2)  the  collision  broadening  of  the  lines. 

In  our  calculations  reliable  values  of  the  line  positions  were  used,  but  the 
collision  broadening  was  only  estimated.  If  the  true  half  width  is  much 
larger  than  the  value  used  here,  absorption  is  enhanced;  if  it  is  much 
smaller,  absorption  is  hindered.  More  accurate  predictions  of  molecular 
absorption  require  reliable  values  for  the  collision  half  widths. 


T [K] 

Fig.  3.  7 Contribution  of  Molecular  Absorption  to  Effective  Absorption 
Coefficient  of  10,  6 Radiation  for  AlgO^  Vapor  in  Equilibrium 
at  One  Atmosphere.  The  coefficient  is  shown  for  three  separate 
COg  laser  lines. 
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Fig.  3.8  Total  Effective  Absorption  Coefficient  for  P(20)  Laser  Line  by  i 

Al^Oj  Vapor  in  Chemical  Equilibrium  at  One  Atmosphere  Pressure. 
The  molecular  absorption  contribution  and  inverse  Bremsstrahlung 
contribution  are  also  plotted. 
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3.  5 Thermal  Conduction  Losses 


The  dominant  loss  at  low  temperature  is  thermal  conduction  of 
energy  from  the  hot  Al^O^  vapor  into  the  cold  surrounding  air.  The 
losses  at  any  given  time  depend  upon  the  local  temperature  profile,  which 
itself  is  determined  by  the  past  history  of  vapor  production,  laser  heating 
and  thermal  conduction  losses.  In  this  report,  however,  a simplistic 
model  has  been  adopted  wherein  there  is  no  significant  heat  addition  or  loss 
during  the  production  of  the  vapor  cloud,  and,  once  the  cloud  is  formed, 
it  heats  without  further  production  of  vapor  from  the  particle. 

The  magnitude  of  the  losses  from  the  vapor  cloud  due  to  thermal 
conduction  are  estimated  by  determining  the  losses  from  a sphere  of 
radius  R(T)  which  is  maintained  at  a constant  temperature  T.  The  losses 
are  therefore  characterized  by  the  thermal  behavior  of  the  air;  that  is, 
by  solving  the  equation 


(p 


c ) 
P 


dT 

dt 


-L  J-  / 

2 Sr  \ 


r2  X 


(3.  15) 


in  the  region  r 3t  R,  where  p c and  X are  the  density,  specific  heat,  and 

pc 

equilibrium  thermal  conductivity  of  air.  In  making  this  approximation,  it 
has  been  assumed  that  the  transport  of  energy  caused  by  the  diffusion  of 
Al^O^  vaPor  into  t^ie  a*r  ant*  t^le  simultaneous  diffusion  of  air  into  the  Al^O^ 
vapor  is  well  approximated  by  the  self  diffusion  within  air  itself.  This 
approximation  seems  reasonable  since  the  diffusion  velocities  of  all  heavy 
components  are  expected  to  be  of  the  same  magnitude. 

The  time  dependent  solution  of  Eq.  (3.  15)  for  constant  Xc  and 
p for  initial  air  temperature  Tq  and  for  vapor  temperature  T^  is  given 
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(T  - T ) = (T  - T ) 
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■(J.  du 


(3.16) 


Ir-Rl 
Zj  Kt 


where  K is  the  diffusivity  \q/q  c^.  The  energy  flux  at  the  surface  of  the 
vapor  cloud  is  simply 


r = R 


X<Ty-To>  ,,  _R_ 

R JilKt 


(3.17) 


For  short  times  the  transient  term  RA/nKt  can  dominate,  but  for  long  times 
the  conduction  losses  approach  a steady  state  value  of  X (Tv-To)/R. 

The  steady  state  limit  can  also  be  determined  for  variable  X c (T). 

The  requirement  that  dT/dt  be  zero  reduces  Eq.  (3.  15)  to 


2 . dT  . . 

r X(T)  ~T~  = constant, 
dr 


(3,  18) 


Invoking  the  conditions  T (R)  = Tv,  and  T(r  ->  ®)  = Tq,  we  find 


X(T)  dT  = R / X(T)  dT 


(3.  19) 


The  flux  at  the  surface  is  then  given  by 

T 

r v 


X(T)  dT 


X dT 
dr 


r = R 


(3.  20) 
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The  quantity  J* XdT  for  air  is  plotted  in  Fig.  3.  9.  The  values  of  X (T)  were 

taken  from  Yos  for  T * 2000°K  and  from  Ref.  3 2 for  2000°K  a T St  300°K. 

Thermal  conduction  losses  are  calculated  as  surface  fluxes.  It  is 

necessary,  however,  to  know  how  those  losses  are  distributed  within  the 

volume  of  the  hot  vapor  cloud.  For  short  times,  the  losses  occur  entirely 

within  a small  layer  of  vapor  adjacent  to  the  cold  air.  At  long  times,  a 

situation  is  reached  whereir.  the  losses  are  spread  uniformly  over  the 

vapor  volume.  Therefore  it  is  important  to  determine  when  the  losses 

affect  the  vapor  cloud  as  a whole,  ratne.  ure  outer  layers.  A 

crude  estimate  has  been  made  as  follows:  The  solution  for  the  cooling  of 

a sphere  initially  at  temperature  T^,  in  the  absence  of  any  hert  sources, 

is  given  in  Ref.  33.  Inspection  of  a plot  of  the  solution  suggests  that  the 

central  temperature  drops  5%  when ./Kt/R2  - . 25.  Assuming  thai  thermal 

-3 

conductivity  of  AhO,  vapor  is  the  same  as  air,  X = 5.  85  x 10  W/cm  K, 

“ ^ o 3°  3 

and  using  p c^  of  the  vapor  at  4000  K (1.3  x 10  J/cm  ),  we  find  the  time 

at  which  the  losses  affect  the  whole  volume  to  be  17pisec  for  the  cloud 

formed  by  evaporating  65.  8%  of  a A1203  particle  of  15 p radius.  This 

time  will  be  used  in  subsequent  calculations  to  determine  when  thermal 

conduction  effects  are  important. 

3.  6 Radiation  Losses 

The  radiation  losses  were  computed  for  a sphere  of  Al-O,  vapor 

u J 

assuming  that  atomic  line  radiation,  atomic  recombination  radiation  and 
Bremsstrahlung  were  the  important  contributions.  This  assumption  is 
valid  at  high  temperatures  where  the  radiation  losses  from  the  plasma 
and  radiative  transfer  to  the  surrounding  air  are  potentially  important 
factors. 

The  absorption  coefficient  for  inverse  Bremsstrahlung  was  deter- 
mined from  Kramer’s  formula  (Eq.  (3.  3))  with  a constant  Gaunt  factor 
of  1.6.  The  calculation  of  the  absorption  coefficient  for  bound-free 

16 

transitions  included  the  transitions  from  all  the  levels  listed  by  Moore. 
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Fig.  3.  9 Integral  of  Equilibrium  Thermal  Conductivity  of  Air  from  300  K to  T 
Steady  s'  ate  conduction  loss  per  unit  volume  from  sphere  of  radius  R 
is  given  by  3 J*  X dT/R2. 
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Estimates  of  the  cross  section  of  photoionization  of  the  ground  state  were 

34 

taken  from  an  survey  of  experimental  data  by  Hudson  and  Kieffer  for 

35 

A1  and  from  the  theoretical  calculation  by  Thomas  and  Helliwell  for  O. 

The  cross  sections  of  all  higher  states  were  determined  by  the  hydrogenic 
36 

approximation.  The  inverse  Bremsstrahlung  absorption  coefficient 

and  the  photoionization  absorption  coefficient  were  added  together  to 

yield  an  effective  continuum  absorption  coefficient.  The  total  surface 

flux  for  continuum  radiation,  including  the  effects  of  reabsorption,  was 

then  determined  for  a sphere  of  uniform  temperature  T and  radius  R(T) 

The  flux  attributable  to  atomic  line  transitions  was  calculated  for 

+ 37  38 

the  lines  and  f numbers  listed  by  the  NBS  tables  for  A1  , A1  and  O.  The 

39 

width  was  determined  by  the  Stark  broadening  parameters  of  Griem.  The 
flux  from  each  line  was  computed  without  including  the  effects  of  other 
lines  or  the  continuum,  but  reabsorption  b"  the  line  itself  was  properly 
included.  The  neglect  of  A1  and  0 atoms*  as  well  a s the  omission  of 
other  broadening  mechanisms, can  be  justified  by  simple  estimates  of 
their  contribution  to  the  overall  flux  in  the  high  temperature  region. 

Stark  broadening  dominates  except  at  low  temperatures  where  radiation 

^ 'f-f'  f 

is  unimportant.  The  O and  A1  ions  may  be  plentiful,  but  their 
emission  occurs  in  the  extreme  ultraviolet  which  is  an  unimportant  region 
of  the  radiation  spectrum. 

The  predictions  for  P ~ . 125,  . 658  and  I.  0 are  plotted  in  Fig.  3.  10. 
The  continuum  radiation  tends  to  be  small  and  transparent.  The  dominant 
losses  are  from  "black”  lines.  The  losses  are  seen  to  be  very  small  at 
low  temperatures;  however,  they  increase  rapidly  and  dominate  over 
thermal  conduction  above  140C0°K. 

The  calculations  discussed  above  are  suitable  for  high  temperatures. 
At  low  temperatures  molecular  emission  may  be  significant.  Fortunately  the 
only  species  which  radiates  in  the  important  region  of  the  spectral  range  at 
low  temperature  is  AlO.  Furthermore  it  is  the  most  prevalent  molecule. 


0 = l.o 
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$ = .125 


w T [10  K] 

Fig.  3.  10  Radiation  Flux  at  Surface  of  Vapor  Cloud  Formed  by  Evaporating 
a Fraction  0 of  a Ai^Oj  Particulate  Which  has  a Radius  of  15 (i. 
Atomic  lines  transitions  and  photo  recombination  are  the  main 
contributions. 
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Therefore  the  contribution  from  all  molecular  emissions  is  expected  to  be 

2 + 2 + 

composed  mostly  of  the  transition  from  the  B £ state  to  the  X £ state 

40 

of  AlO.  The  oscillator  strength  for  all  transitions  (Av  = 0,  _+  1)  is 

- 2 - 1 

1. 12  x 10  . The  band  occurs  near  uu  = 20635  cm  . Assuming  trans- 

parent emission, which  is  the  maximum  possible  emission,  we  find  the 
radiation  loss  from  molecular  radiation  to  be  210  W/cm  at  4000°K.  This 

6 2 

loss  is  less  than  two  percent  of  the  absorption  rate  for  I = 3 x 10  W/cm  , 
thus  molecular  emission  is  unimportant  for  the  determination  of  heating 
rates. 

3.  7 No  Loss  Approximation 

Shortly  after  vaporization  begins,  a layer  of  vapor  is  produced 
which  protects  the  vapor  subsequently  produced  from  the  cooling  effect 
of  thermal  conduction  into  the  air.  In  Subsection  3.  5 the  time  for  which 
this  prelection  lasts  was  crudely  estimated  to  be  17|jsec  for  a 65.  8% 
vaporized  particle.  Therefore  the  first  few  microseconds  of  vapor  heating 
occur  uninfluenced  by  losses  due  to  thermal  conduction.  The  "no  loss" 
approximation,  which  is  valid  at  short  times,  is  described  by 


dT  ' 

p Cp  * = KL  - 


(3.21) 


and  the  time  to  heat  the  vapor  to  temperature  T starting  from  an  initial 


temperature  Tg  is  found  from 


-tf 


P 

-T2-  dT 


K 


(3.22) 


A plot  of  T versus  I t for  T = 4000°K  is  given  in  Fig.  3.  11.  Since 
r o s 

similar  plots  for  different  choices  of  T correspond  to  shifts  in  the  origin 

s 

of  the  fluence  co-ordinate*  negative  values  of  fluence  were  included.  The 


; < 


fj 


■■  ^ ..  | 


-55- 


'Sts?****** 


Io  t [j/cm2] 

Fig.  3.  11  Temperature  of  Equilibrium  Vapor  as  a Function  of  Fluence  - 

No  Loss  Approximation.  "X"  is  used  to  designate  each  one  hundred 
degree  temperature  point  from  3900  K to  4500  K. 
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graph  dramatically  demonstrates  that  most  of  the  time  spent  heating  the 

vapor  occurs  during  the  initial  heating  to  5000°K.  In  the  "no  loss" 

approximation,  the  fluence  required  to  heat  the  vapor  from  4000°K  to 

20000  K is  approximately  71  J/cm  . For  comparison  recall  that  the 

2 

fluence  necessary  to  initiate  vaporization  is  41.4  J/cm  while  the  fluence 

needed  to  vaporize  a fraction  3 of  the  original  particle  having  a 15|i 

2 

radius  is  225  3 J/cm  . 

If  the  heating  always  occurred  without  losses,  the  time  to  vaporize 

a particle  and  then  to  heat  the  cloud  to  20000°K  would  scale  as  IQt.  This 

scaling  breaks  down  at  low  intensity  for  two  reasons.  The  first  reason, 

discussed  in  Section  2,  is  that  a pure  vapor  cloud  cannot  be  produced. 

A second  reason  is  that  the  vapor  does  not  heat  sufficiently  rapidly  to  allow 

o 

louses  to  be  neglected.  The  time  to  heat  the  vapor  to  20000  K ’ s very 

sensitive  to  the  temperature  the  vapor  cloud  has  achieved  when  losses 

become  important.  Using  17psec  as  the  time  over  which  the  "no  loss" 

approximation  is  valid,  we  see  from  Fig.  3.  11  that  an  intensity  of  3 x 10^ 

2 o 

W/cm  heats  the  vapor  to  4300  K before  losses  must  be  included.  If  the 

5 2 o 

intensity  is  5 x 10  W/cm  , however,  the  vapor  reaches  20000°K  before  the 

17|isec  has  passed.  Thus  it  is  only  at  the  lowest  intensities  studied  here 

that  thermal  conduction  losses  become  important. 

It  should  be  noted  that  the  model  of  constant  pressure  heating  in 

thermodynamic  equilibrium  does  not  accurately  describe  the  rapid  heating 

6 2 

at  high  temperatures.  Even  for  a flux  of  3 x 10  W/cm  the  model  predicts 
that  the  vapor  heats  from  10000°K  to  20000°K  within  . 3jaseconds  - a time 
of  the  same  magnitude  as  the  time  required  to  equilibrate  pressure  and 
maintain  chemical  equilibrium.  However,  since  the  heating  rate  at  high 
temperature  must  be  slowed  down  by  a factor  of  ten  in  order  alter  the 
total  heating  time  significantly,  the  possibility  of  a non -equilibrium  vapor 
is  not  pursued  here. 
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3.  8 Averaged  Los  a Approximation 

Although  the  no  loss  calculations  are  applicable  for  strong  beams,  such 

6 2 

as  8 x 10  W/cm  , as  the  intensity  is  lowered  the  time  to  heat  the  vapor 
becomes  longer  than  the  maximum  time  over  which  conduction  can  be  neglected. 
A simple  method  of  accounting  for  losses  is  to  average  them  over  the  entire 
volume  of  the  vapor.  In  earlier  subsections  J*  X dT  and  have  been 

determined.  The  temperature  history  of  the  vapor  cloud  in  an  averaged  loss 
approximation  is  described  by 


P c 


dT  1 

= K I 

La  O 


p dt 


3 r X dT  3 SRAD 

R2  " R 


(3.23) 


The  first  calculation  to  be  made  in  the  model  is  the  threshold  intensity  I 

1 ri 

below  which  the  laser  cannot  heat  the  vapor.  The  right  hand  side  of  Eq.  (3.23  ) 
must  be  positive  for  heating  to  occur,  therefore  the  quantity 


I * 

s 


K R 
L 


+ s 


RAD 


is  the  laser  Intensity  which  would  just  sustain  the  vapor  cloud  at  a given 

temperature.  A plot  of  I permits  the  threshold  intensity  for  heating  of 

the  vapor  cloud  to  be  determined  by  inspection.  This  quantity  is  plotted 

in  Fig.  3.  12  for  three  choices  of  8> 

In  no  heating  were  to  occur  without  losses  being  important,  the 

6 2 

threshold  intensity  is  1,9  x 10  W/cm  for  an  entirely  vaporized  particle. 

-2/3 

The  threshold  scales  as  8 for  partially  vaporized  particles.  Of  course 
it  takes  infinite  time  to  achieve  breakdown  at  threshold.  Experiments 
limited  to  a 100|isec  pulse  length  should  observe  higher  thresholds,  since 
the  breakdown  must  occur  while  the  laser  is  on. 
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In  order  to  compare  the  trends  predicted  in  this  section  with 

experimental  observations,  we  focus  on  vapor  clouds  with  9=65.  8% 

since  these  are  the  largest  clouds  expected  to  be  produced  within 

100 Usee.  From  Fig.  3.  12  the  threshold  intensity  is  predicted  to  be 

4.  3 x 10  W/cm  if  losses  are  always  important.  In  fact,  however, 

some  no  loss  heating  does  occur,  and  the  threshold  is  lowered.  At 
6 2 

4 x 10  W/cm  we  expect  no  loss  heating  to  occur  for  ~ 17  ^isec  at  which 

time  the  typical  temperature  will  be  5000 °K.  There  is  then  no  problem 

heating  the  cloud  from  this  temperature  to  20000°K,  even  when  losses  are 

6 2 

included,  A more  reasonable  estimate  of  the  threshold  is  3 x 10  W/cm6, 

No  loss  heating  raises  the  temperature  to  4300°K.  Then  the  vapor  can 
heat  at  a rate  determined  by  the  averaged  loss  model.  The  additional 
time  to  heat  to  temperature  T is  determined  by 

T 

-/ 

4300 

6 2 

This  number  has  been  calculated  for  P = 65.8%  and  I = 3 x 10  W/cm  . 

o 

The  primary  difference  between  the  no  loss  time  and  the  averaged  loss 
time  occurs  during  the  heating  to  5000°K.  The  heating  takes  5 ^sec  in 
the  no  loss  calculation,  but  it  requires  27  Usee  in  the  averaged  loss  model. 

The  total  time  to  heat  the  cloud  is  17  ^sec  to  reach  4300°K,  another  27  ^sec 
to  reach  5000°K  and  2^  sec  for  all  subsequent  heating.  The  total  time,  including 
the  time  to  produce  the  vapor,  is  60  ^sec  if  heating  of  the  vapor  commences 
as  soon  as  vaporization  begins,  but  it  is  96  iisec  if  the  heating  does  not  occur 
until  the  full  65.  8%  of  the  particle  is  vaporized. 
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Fig.  3.  12  Plot  of  Laser  Intensity  Required  to  Offset  Losses  for  Three 
Choices  of  8,  the  Percent  of  a 30,  p Diameter  Particle 
which  is  Vaporized. 
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In  order  to  demonstrate  the  Importance  of  losses  upon  the  heating 

rate,  we  have  plotted  in  Fig.  3. 13  the  time  rate  of  change  of  temperature 

for  two  laser  intensities.  The  heating  rates  for  the  "no  loss"  approximation, 

as  well  as  for  the  "averaged  loss"  approximation  with  3 = 1.  0 and  3 = • 658, 

6 2 

are  shown.  The  losses  have  only  a small  effect  for  I = 8 x 10  W/cm  , 

6,2° 

but  they  are  very  important  for  Iq  = 3 x 10  W/cm  . Indeed,  the 

temperature  at  which  the  transition  is  made  from  "no  loss"  to  "averaged 

loss"  is  crucial  in  determining  whether  or  not  breakdown  occurs  for  3 = * 658 

and  I = 3 x 106  W/cm2. 
o 

In  summary  the  models  for  vapor  heating  predict  that  the  vapor  heats 

6 2 

in  times  consistent  with  the  observed  delay  times  for  I a 5 x 10  W/cm  . 

6 °2 

The  models  also  predict  a threshold  near  3 x 10  W/cm  , but  the  delay 
time  near  the  threshold  cannot  be  predicted  unambiguously. 

If  the  laser  wavelength  were  changed  to  3.8|i,  molecular  absorption 
would  not  be  significant  and  inverse  Bremsstrahlung  absorption  would  drop 
by  a factor  often.  The  threshold  intensity  would  therefore  increase  by  at 
least  a factor  often.  However,  it.  is  possible  that  other  particulates  may 
produce  vapors  which  absorb  readily  at  3.  8ji,  and  which  therefore  may 
heat  rapidly  at  intensities  well  below  the  expected  threshold  for  A1203. 

Each  particulate  and  each  laser  wavelength  must  be  studied  individually. 


Fig.  3. 13  Rate  of  Change  of  Temperature  with  Time,  as  a 


Function  of  Temperature  for  Equilibrium  AI2O3 
Vapor.  Results  for  three  models  of  losses  and  two 
laser  intensities  are  shown. 
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4.  CREATION  OF  AN  AIR  PLASMA 


A single  particulate  cannot  cause  breakdown  without  first 

creating  an  air  plasma,  since  the  vapor  remains  transparent  to  laser 

o * 

radiation.  Recall  that  at  20000  K the  absorption  coefficient  K is 

-1  L 
approximately  . 6 cm  and  the  radius  R for  an  entirely  vaporized  particle 

is  only  . 13  cm.  Thus  the  maximum  optical  depth  is  , 16.  In  this  section 

a simple  estimate  of  the  time  to  create  an  air  plasma  is  made. 

First,the  temperature  at  which  the  air  efficiently  absorbs  the  laser 

energy  must  be  determined.  We  adopt  the  terminology  used  in  laser 

supported  combustion  wave  models  and  call  this  temperature  the  ignition 

temperature  T^.  For  air  the  only  absorption  mechanism  at  10.  6 n is 

inverse  Bremsst  rah  lung.  In  Fig.  4.  1 a plot  of  the  absorption  coefficient 

of  air  is  compared  to  the  absorption  coefficient  for  Al2C>3  vapor.  The 

air  absorption  rate  equals  the  vapor  absorption  rate  at  14,  000°K.  When 

the  air  reaches  the  14000°K,  therefore,  it  is  as  effective  as  the  vapor  in 

attenuating  the  laser  beam.  At  intensities  of  interest,  however,  the  air 

may  heat  rapid*/  at  considerably  lower  temperatures.  A plot  of  dT/d(tIQ) 

for  air  is  shown  in  Fig.  4.  2.  Rapid  heating  occurs  when  the  rate  of  change 

of  temperature  is  10000°K/iasec  = 10^  K/ sec.  For  a laser  where  intensities 
6 2 

are  near  5 x 10  W/cm  the  value  of  dT/d(t  I ) should  therefore  be  of  the 
3 2 ° 

order  of  2 x 10  Kcm  /J.  On  the  basis  of  the  heating  rate,  the  ignition 
temperature  may  he  as  low  as  12000°K. 

For  simplicity  we  assume  that  the  Al^O^  vapor  is  at  20000°K.  It 
transfers  energy  to  the  surrounding  air  and  thereby  heats  a layer  of  air 
to  l,.  Once  the  temperature  of  the  air  reaches  T.,  it  is  rapidly  heated  by 
the  laser  to  a high  temperature  - assumed  to  be  20000°K  in  the  calculations 
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Fig.  4.  1 Effective  Absorption  Coefficiei 


Vapor  and  for  Air. 
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Fig.  4.  2 The  Rate  of  Change  of  Temperature  with  Finance  as  a Function 
of  Temperature  for  10.  6 H Radiation. 
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performed  in  this  election  - after  which  the  newly  heated  air  can  in  turn 
transfer  energy  to  the  surrounding  layers,  of  cold  air.  The  rate  of 
expansion  of  the  plasma  is  given  approximately  by 


pahay^q 


(4.  1) 


where  p is  the  density  of  air  at  one  atmosphere  and  20000°K,  Ah  (T^)  is 

the  energy  required  to  heat  air  from  300°K  to  and  q is  the  energy  flux 

from  the  hot  plasma.  An  assumption  implicit  in  Eq.  (4.  1)  is  that  all  the 

flux  q is  effective  in  heating  the  air  in  the  layer  adjacent  to  the  plasma. 

The  flux  q has  two  components.  For  small  radii,  less  than  . 05  cm, 

4J 

thermal  conduction  is  known  to  be  dominant.  For  large  radii  the  energy 

42  43 

flux  is  controlled  by  radiation.  ’ Hera  we  are  dealing  with  a transition 
zone  in  which  both  mechanisms  contribute.  In  order  to  demonstrate  that 
the  air  is  rapidly  heated  we  specialize  to  a specific  example.  We  choose 
a particle  which  is  65.  8%  vaporized.  From  Fig.  3.  5 the  vapor  cloud 

is  . 09  cm  in  radius  at  the  start  of  the  air  heating.  The  radiative  flux 

2 

in  the  vacuum  ultraviolet  (VUV)  is  approximately  5400  W/cm  . Of  this 

2 

radiation  only  200  W/cm  lies  in  the  range  hv  i 11  eV  which  is  the,  spectral 
region  absorbed  by  oxygen  and  nitrogen  atoms.  The  rest  of  the  radiative 
flux  lies  in  the  Schumann- Runge  band  of  molecular  oxygen.  The  fLux  due 
to  thermal  conductivity  is  given  by  (see  Eq.  (3.  17)) 


q.AS 

q R 


Vn  Kt 


which,  for  air  at  20000  K and  t = 10  sec,  is  1.  5 x 10  /cm  . Therefore 

4 2 

the  total  flux  into  the  air  from  the  hot  Al^C^  vapor  is  2.  0 x 10  W/cm  . 
On  the  other  extreme,  once  the  plasma  reaches  . 5 cm  in  radius,  the 
radiative  flux  in  the  region  hv  k 11  eV  is  enhanced  and  becomes  nearly 
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3.  8 x 10  W/chi  so  that  the  total  radiation  ia  4.  3 x 104  W/cm  . The  flux 

due  to  thermal  conduction  is  small  because  of  the  preheating  caused  by- 

radiation.  As  a first  approximation,  therefore,  a value  of  q = 2,  0 x 10^ 

2 

W/cm  is  Ufod  to  Characterize  the  full  regime  of  interest.  This  leads  to 
a slight  overestimate  of  the  time  necessary  to  produce  large  plasmas. 

The  velocity  of  the  plasma  front  can  now  be  estimated.  For  of 
14000°K,  we  find  dR/dt  = 5 x 10^  cm/soc*  and  for  of,  12000°K,  (^R/dt 
increases  to  7 x 10^  cm/sec.  At  either  rate  the  plasma  doubles  in  radius 
in  less  than  two  microseconds.  If  these  estimates  of  expansion  remain 
accurate  for  radii  as  large  as  .5  cm,  it  takes  a maximum  of  8 microseconds 
to  produce  a plasma  1 cm  in  diameter  (and  perhaps  the  time  is  as  low  as 
6(isec).  The  purpose  of  these  calculations  is  to  demonstrate  that  the  air 
plasma  forms  very  rapialy  after  the  vapor  reaches  high  temperatures. 

The  details  of  the  growth  of  an  air  plasma  are  discussed  in  the  next  section. 

The  effect  of  laser  wavelength  and  particle  composition  upon  this 
stage  of  the  breakdown  process  is  uncertain.  The  air  heats  at  a slower 
rate  if  the  wavelength  is  shortened.  However,  at  the  intensities  considered 
here,  the  heating  is  still  rapid.  Since  the  laser  is  not  directly  responsible 
fer  the  heating  of  th<-.  cold  air,  the  most  important  parameters  are  the 
temperature  and  size  of  the  plasma  core,  not  the  laser  intensity.  Of  course, 
the  temperature  of  the  core  is  related  to  the  laser  intensity,  but  for  the 
intensities  required  to  initiate  vapor  heating  the  core  temperature  of  the 
resultant  plasma  is  expected  to  exceed  20, 000°K.  In  short,  changes  of 
wavelength  and  particle  composition  are  not  expected  to  alter  the  estimates 
of  this  section. 
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5.  AIR  PLASMA  DYNAMICS 


The  previous  sections  have  dealt  with  the  absorption  of  laser 
energy  by  solid  particles,  which  are  either  completely  or  partially 
vaporized,  and  the  transfer  of  the  energy  from  the  vapor  to  the  surround- 
ing air.  The  hot  air  is  ultimately  heated  via  inverse  Bremsstrahlung 
absorption  until  the  expanding  air  plasma  becomes  opaque  to  the  laser 
radiation.  This  p’asma  propagation  mechanism  is  similar  to  that  of  the 
laser- supported  combustion  (LSC)  wave  propagation.  For  a given 
particulate  distribution,  it  is  difficult  to  estimate  the  corresponding  initial 
air  plasma  bubble  distribution  for  a given  pulse  time  and  laser  intensity 
because  the  particulates  may  not  be  completely  vaporized.  Hence,  in 
this  section,  we  will  examine  the  air  plasma  growth  dynamics  with  the 
initial  size  assumed  to  be  given.  For  simplicity,  we  will  also  ignore  the 
presence  of  particle  vapor  which  may  be  important  in  some  cases. 

Initial  Conditions 

The  late-time  air  plasma  growth  begins  when  sufficient  energy  has 
been  transfered  from  the  vapor  to  the  surrounding  air  such  that  the 
high  temperature  air  plasma  begins  to  absorb  the  laser  energy  directly. 

The  initial  size  of  the  air  plasma  depends  obviously  on  the  size  of  the 
vapor  which,  in  turn,  depends  on  the  original  particulate  size,  the  fraction 
of  solid  that  vaporizes  and  the  laser  intensity. 

As  mentioned  before,  the  particulates  of  interest  are  primarily 

A1„0_,  SiO  , and  graphite.  Now,  consider  the  diameter  of  the  solid  parti- 
es 3 c. 

culates  to  be  10  - 20  jjm  and  the  pressure  to  be  1 atm.  Assuming  that  the 
solid  particulates  vaporizes  at  1 atm  and  vaporization  temperature  and 
that  the  equilibrium  composition  prevails  at  T^,  one  obtains  the  vapor  size 
(for  complete  vaporization)  to  be  several  hundred  (im,  Table  5.  1. 


Table  §.  1 


Particulates 


A12°3 


Graphite 


Tv*  K 
3200 
2500 
5100 


Components  at  T Vapor  Diameter,  ^m 


O,  Al,  AlO 
SiO,  O, 


30C  - 600 
230  - 460 
420  - 840 


The  particulates  may  not  be  completely  vaporised  and  the  vapor  diameters 

would  then  be  less  than  the  above  values.  Assuming  complete  vaporisation, 

we  arbitrarily  take  the  initial  air  plasma  sise  to  be  twice  the  average  vapor, 

i.  e. , the  initial  radius  = 0.  05  cm.  As  will  be  shown  later,  the  initial  radius 

will  be  parameterised.  Finally  from  the  effective  laser  absorption  coeffi- 

43 

cient  for  10.  6}j  radiation  in  air  at  1 atm,  one  finds  that  the  absorption 

-2-1  -1 
coefficient  rises  from  2x10  cm  at  10,  000  K to  approximately  0.  8 cm 

at  15000  K.  For  simplicity,  we  choose  the  initial  temperature,  T^  = 14,000, 

which  can  also  be  parameterised  to  assess  the  effects  of  T^  on  the  final 

solution. 

Laser  Absorption  and  Plasma  Emission 

The  phenomenology  of  air  plasma  growth  involves  coupled  fluid 

dynamics  and  radiant  energy  transfer.  As  the  air  plasma  grows,  the 

energy  transfer  mechanism,  which  is  analogous  to  that  of  laser  supported 

combustion  (LSC)  waves,  is  dominated  by  radiation.  Since  radiation  is 

expected  to  control  the  maximum  temperature  and  the  expansion  velocity 

of  the  air  plasma,  it  is  important  to  properly  characterize  the  radiative 

transport.  The  local  radiative  properties  are  determined  by  the  effective 

43 

absorption  coefficient  for  inverse  Bremsst  rah  lung 


Kv  = * 


16tl2  ( JJL  -\ 

3 \3  mek  ) 


1/2  ,26  (l  - e"“/kT)n 

5 — rrr — [n  + 4 n ] 

c m ,3  _l/2  1 + ++J 


(U)L)  T 


(5.  1) 


and 

g = 0.55  in  27(^4/3(pxe)"1/3 

where  k is  Boltzmann's  constant; -ft  is  Planck's  constant;  c is  the  speed  of 

light;  m is  the  electron  mass;  ta.  is  the  laser  photon  energy;  T is  the 

temperature;  p is  the  pressure;  X is  the  mole  fraction  of  electron;  and 

n , n. , n,  , are  the  number  densities  of  electrons,  ions  and  double  ions 
• + ++  43 

respectively.  Pirri  et  al,  have  shown  that  the  effective  absorption 
coefficient  at  1 atm  is  approximately  0.  8 for  temperatures  greater  than 
14000  K (cf  Fig.  5. 1).  The  initial  diameter  of  the  air  plasma  is  of  the 
order  of  0. 1 cm  or  less.  Hence  the  source  term  in  the  energy  equation 
can  be  assumed  to  be  volumetric  absorption  i.  e. , K^I  where  I is  the 
laser  intensity. 

43 

For  the  radiative  emission  from  the  hot  air  plasma,  Pirri  et  al, 
have  shown  that  it  can  be  divided  into  the  transparent  (iu  < 11  eV)  and  the 
black  (u)>  11  eV)  contributions.  The  transparent  band  can  be  written  as 


= 6.  3 x 10”10  T n 

o 


kT 

I . 
air 


+ 4.  2 


-I  . /kT 
air 
e 


(5.2) 


where  n is  the  number  density  of  neutrals;  1 . is  the  ionization  potential 
o air 

for  air  and  q+/qQ  is  the  ratio  of  the  partition  functions  of  ions  to  neutrals. 
The  UV  band  (iu  >11  eV),which  accounts  for  most  of  the  axial  transport  is 
essentially  black,and  given  by 


where 


I f 

II  eV  11 eV 


U) 


u> 


kT 


(5.  3) 
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is  the  black  body  function  and  v is  the  frequency.  For  the  present 
calculations,  we  find  that  the  cor* ribution  from  the  transparent  band  is 
negligibly  small  compared  with  the  UV  band  and  hence  the  plasma 
emission  in  the  following  analysiu  will  be  represented  by  the  UV  band 
alone.  With  the  radiative  transport  model  adapted  as  above,  we  can  now 
proceed  to  examine  the  air  plasma  dynamics  by  means  of  a perturbation 
analysis. 

Perturbation  Analysis 

Consider  a spherical  air  plasma  of  radius  R.  The  governing  equation 
for  the  growth  dynamics  can  be  written  as 


+ ilcui  + 2p_u  _ Q 
bt  br  r 

(5.  4) 

bu  _bu  _ 2. 

bt  U br  p br 

(5.  5) 

bs  . bs  V 3tt 
0-sr  • pub?=~-T5: 

(I  + ar) 

(5.6) 

p = (y  - 1)  pY  exp  I ) 

(5.7) 

V v / 

2 

v ' C x3/2 

(5.8) 

dR  2nkv3T  -hv/kT 
dt  " ...  2 * 

(5.9) 

where  p is  tho  density;  u is  the  radial  velocity;  S is  the  entropy;  y is  the 
ratio  of  specific  heats;  is  the  heat  capacity  at  constant  volume;  Ah%  is 
the  enthalpy  difference  between  the  hot  air  at  the  edge  of  air  plasma  and 
the  cold  ambient  air;  the  subscript  a indicates  ambient  condition;  and  t 
and  r are  time  and  radial  distance  respectively.  Equation  (5. 9)  is  derived 
by  equating  the  UV  band  emission  to  the  energy  consumed  by  the  ingestion 
of  cold  ambient  air  which  is  then  brought  to  the  air  pUtsma  temperature. 
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The  above  Eqs.  (5.4)  - (5.9)  can  be  normalised  by  the  following  characteristic 


quantities;  the  initial  sound  speed,  a - */yRT~,  the  initial  air  plasma  radius, 
Rj  and  a heating  time,  T.  Then  the  nondimensional  variables  are 


— p — u “ t — 

* * fr;  u=i 1 e j p 


and  the  conservation  equations  become 


£o-+  g [±iLsl  + i^-1  = o 

^ L dr  r J 

f^  + s ; %1  = » 

* l »r  P Jrj 

\6ttv  kT  -uj./T 

PT  **-  + e Bi»  ~r  = «IS,-(1+A-) — ~ae  <5- 

SF  dr  \ a a/  1 C R 


(5.  10) 


(5.11) 


12) 


«£_  -v5 

* ' ^ pm  C2  . 


(5.  13) 


The  characteristic  heating  time  T has  been  defined  as  T ■ T^  Rj/I. 

Two  basic  parameters  appear  in  these  normalised  equations,  the  Bouguer 

Number,  a * R^  which  is  the  ratio  ot  the  air  plasma  radius  to  the 

absorption  length  at  the  laser  frequency,  and  the  Boltsmann  Number, 

S = o.  C T„  a/I,  which  is  the  ratio  of  the  convective  energy  flux  to  the 
Ki  v i 

incident  radiation  energy  flux. 

6 7 2 

For  the  laser  intensity  of  10  - 10  W/cm  , the  initial  temperature 
of  14000  K and  the  initial  pressure  of  1 atm,  the  Boltsmann  Number  becomes 
much  less  than  unity.  If  each  of  the  dependent  variables  in  the  nondimen- 
sionalised  equations  is  expanded  in  a series  of  functions  whose  coefficients 
are  powers  of  B,  i.  e. , 
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— 2 
p * pQ  + 0 Pj  + fi  p2  + 

the  equations  can  be  reduced  to  sets  of  equations  which  appear  as  coefficients 
of  0.  The  sero-order  solution  of  the  nondimensional  momentum  equation 
provides 

u = 0 
o 

which  indicates  that  there  is  no  flow  and  that  the  growth  of  the  air  plasma 
is  by  LSC  wave  mechanism  alone.  As  a result,  the  problem  is  reduced  to  a 
solution  of  the  energy  equation  and  the  LSC  wave  propagation  equation. 


bT  „ ,, 

ot  o 


(t  t-a_\ 

\ 4ha  / I C2 


dR  2nv3kp.  C T2  -uu./T 

o x v i __  i o 

— . T e 

ot  ,,  p 2 _ o 

Ah  p C I 
a Ha 


(5.  U) 


(5.  15) 


Equations  (5.  14)  and  (5.  15)  have  been  solved  numerically  fox*  Tj  = 14,  000  K, 

= 1 atm,  Ri  = . 05  cm  and  laser  intensity  of  10^,  5 x 10^  and  10^  W/cm2. 

With  the  LSC  mechanism  alone,  the  temperature  typically  climbs  from  to 

an  almost  plateau  temperat  ure  in  less  than  one  ^second.  An  examples  tempera- 

6 2 

ture  profile  for  Iq  = 5 x 10  W/cm  is  shown  in  Fig.  5.  2.  The  air  plasma  growth, 
on  the  other  hand,  is  minimal  as  shown  in  Fig.  5.  3,  For  a lOOpigea  pulse,  the 
final  plasma  radii  at  the  end  of  the  pulse  are  only  0.  057,  0.  072  and  0.  085  cm 

c 6 7 2 

for  laser  intensities  of  10°,  5 x 10  and  10  W/cm  respectively.  However, 
as  the  intensity  increases,  the  average  pressure  inside  the  ai  r plasma 
increases  rapidly.  For  example,  at  10  W/cm  it  is  estimated  that  the 
average  over-pressure  reaches  100  atm  at  an  early  time.  At  the  same  time, 
the  characteristic  time  for  pressure  equilibration  at  14, 000  K is  estimated 
to  be  0.  3 Msec  which  is  short  compared  with  the  pulse  time  of  interest,  i.  e. , 
lOOpsec.  Hence,  the  air  plasma  growth  becomes  dominated  by  the  hydro  - 
dynamic  effects.  In  the  following  section,  we  will  develop  a model  which  will 
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incorporate  both  the  ingestion  of  mass  by  the  LSC  wave  mechanism  and  the 
hydrodynamic  expansion  as  a result  of  the  ingestion. 

The  "Snow  Plow"  Analysis 

Consider  a spherical  air  plasma  which  absorbs  laser  energy.  The 

heated  zone  propagates  outward  with  the  same  radiative  transfer  mechanism 

as  the  LSC  wave.  This  outward  propagation  ingests  high  density  ambient 

air.  As  a result  of  the  mass  ingestion,  the  average  pressure  rises. 

Here  it  is  assumed  that  the  pressure  is  instantaneously  equilibrated  and 

hence  the  density  and  the  temperature  are  uniform  inside  the  air  plasma. 

This  high  pressure  air  plasma  is  then  allowed  to  expand  and  "push”  away 

the  ambient  air  in  a manner  similar  to  that  of  the  "snow  plow".  For  simplicity, 

the  velocity  profile  inside  the  air  plasma  is  assumed  to  be  always  a linear 

profile,  i.  e. , u ~ u (t)  r/R  where  u is  the  radial  velocity  at  the  air 

R R 

plasma  edge. 

From  the  conservation  of  mass,  one  gets 

dRLSC  }q_  dR_ 
dt  R dt  ' K dt 


The  first  term  on  the  RHS  represents  the  increase  of  density  due  to  the 
ingestion  of  mass  while  the  second  term  indicates  the  decrease  of  density 
due  to  the  air  plasma  expansion.  Here  the  expansion  equation  includes  the 
flow  velocity  and  the  LSC  wave  contribution. 


dR  _ dRLSC  _ EB 

dt  ~ UR  dt  UR  p Ah 

a a 

Hence  continuity  becomes 


(5.  16) 


d£  _ 2 

dt  R 


(5.  17) 


From  the  momentum  equation 


one  can  integrate  over  the  complete  volume  of  the  air  plasma.  Remembering 
p = p(t)  and  Gauss’s  theorem,  we  obtain 


R(t) 
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bu  2 , i. 

— 4irr  dr  +- 
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4n  r^  dr 


(p-  V 
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4TTR~ 


With  the  linear  velocity  profile  assumption,  the  momentum  equation  can 
finally  be  written  as 


&UR  4 (p  ~ pa) 

bt  = Rp 
For  the  energy  equation,  one  has 


u E_ 
R B 

R p Ah 

a 


where  e ie  the  internal  energy  and  Q represents  the  radiative  transport 
(the  laser  absorption  and  the  plasma  emission).  Again,  we  average 

the  entire  energy  equation  over  the  air  plasma  volume  by  integrating.  Using 
the  continuity,  momentum  and  the  expansion  equations,  one  gets 


be 

bt 


5 p R 


(4  P. 


9 p)  + E& 
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R p Ah 


f 9 2 
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(5. 19) 
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where  E is  the  laser  absorption  term  which  can  readily  be  derived  for 
ft 

uniform  temperature  to  be 


e .-ai- 

a 2 p R 


J. 

2 


2 R K 


•2RKV 


1 
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andthe  internal  energy  can  be  expressed  as 

e = h(T)  - £ 

P 

Finally,  Eqs.  (5,  16)  - (5.  19)  together  with  the  equation  of  state 


p = pRT  (5.  20) 

provide  five  equations  for  the  five  unknowns  p,  T,  P,  R and  u . The  enthalpy 

44  R 

h = h (T,  P)  is  obtained  from  the  results  of  Yob.  The  above  equations 

have  been  solved  numerically  for  an  initial  temperature  of  14,  000  K,  initial 

6 6 7 2 

pressure  of  1 atm  and  laser  intensities  of  10  , 5 x 10  and  10  W/cm  . As 
shown  in  Fig.  5.4  for  the  case  of  R^^  = 0.  05  «m,  the  temperature  distributions 
are  not  significantly  different  between  the  "snow  plow"  model  and  the  calcula- 
tions considering  only  the  LSC  mechanism,  because  the  energy  equation  is 
dominated  by  the  radiative  transports.  However,  the  air  plasma  growth 

rate  has  increased  a great  deal  due  to  the  hydrodynamic  effects.  For 

6 2 

R^  = 0.  05  cm  and  I = 5 x 10  W/cm  , the  air  plasma  radius  doubles  in 
approximately  1 \isec,  (cf  Fig,  5.  5).  One  notes  that  the  acceleration  occurs  in 
an  early  time  and  for  each  intensity  the  velocity  reaches  a different  plateau. 

For  all  practical  purposes,  we  can  use  the  plateau  velocity  for  each  intensity 
to  assess  the  air  plasma  growth  rate  in  the  pulse  time  of  interest. 

The  plateau  values  of  the  expansion  velocity  are  plotted  against  the  laser 
intensity  in  Fig.  5.6,  Since  the  particulates  maybe  only  partically  vaporized, 
calculations  have  also  been  done  for  R^  = 0,  025,  0.  010  and  0.  005  cm.  For 
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Expansion  Velocity,  cm/sec 


comparison,  the  results  for  R^  = 0.  05  cm  from  the  perturbation  calculations 
are  also  presented.  It  can  be  observed  that  although  the  LSC  mechanism 
is  vital  in  ingesting  the  high  density  ambient  air,  the  air  plasma  growth 
is  primarily  a hydrodynamic  effect,  and  that  for  a given  pulse  width  a decrease 
in  R^  requires  an  increase  in  laser  intensity  in  order  to  reach  a given  final 
air  plasma  size.  When  the  expansion  velocity  passes  a value  of  about 

A 

3 x 10  cm/sec, which  is  the  sound  speed  in  the  ambient  air,  the  results 

imply  that  there  will  be  shock  wave  proceeding  the  air  plasma  front.  This 

phenomenon  has  not  been  accounted  for  in  the  present  analysis.  Consider 

an  arbitrarily  chosen  final  air  plasma  size  to  be  1 cm  in  diameter.  Figure  5,  7 

shows  the  intensity  requirement  for  a given  initial  size  air  plasma  to  grow 

to  1 cm  in  diameter  in  50  and  lOO^Jlsec.  Suppose  that  we  have  a 20  cm 

6 2 

diameter  laser  beam  with  an  intensity  of  4 x 10  W/cm  and  the  particulate 
distribution  is  such  to  provide  an  average  intitial  air  plasma  bubble  size  of 
R = . 025  cm.  Then,  when  the  number  density  of  these  bubbles  is  greater 

i _2 

than  1.  3 cm  across  the  beam,  the  air  plasma  will  be  able  to  fill-up  the 
beam's  cross  section  in  50psec. 


Laser  Intensity,  10  W/cm 


Fig.  5.  7 The  Required  Laser  Intensity  for  an  Air  Plasma  to  Reach  1 cm 
in  Diameter  in  a Given  Time  (50  or  lOO^isec)  as  Function  of 
Initial  Radius. 
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6.  SUMMARY 


The  features  of  long  pulse  breakdown  in  particulate  contaminated 
air  which  distinguish  it  from  short  pulse  breakdown  ares  (1)  low  threshold 
intensities  of  the  order  of  several  megawatts  per  square  centimeter, 

(2)  long  delay  times  of  the  order  of  tens  of  microseconds,  and  (3)  the 
production  of  a luminous  plasma  several  microseconds  before  breakdown 
occurs.  These  features  suggest  that  breakdown  occurs  via  heating  of 
the  particulate  vapor  at  constant  pressure  and  in  thermodynamic  equili- 
brium, with  the  subsequent  production  of  an  air  plasma  through  energy 
transfer  from  the  hot  vapor.  This  mechanism  is  studied  theoretically 
in  the  four  preceding  soctions,  and  the  results  of  the  investigation  are 
consistent  with  experimental  observations. 

As  an  example,  consider  an  alumina  particle  having  a diameter  of 

6 2 

30p.  When  this  particle  is  subjected  to  5 x 10  W/cm  of  10-  6 radiation, 
we  predict  the  following  time  scales:  (1)  The  particle  is  50%  vaporized  within 
23  microseconds,  (2)  The  luminous  vapor  thus  formed  heats  without  losses 
and  reaches  20000°K  within  14  microseconds,  (3)  The  hot  vapor  cloud 
transfers  energy  to  the  surrounding  air,  thereby  enabling  the  air  to  absorb 
laser  radiation.  The  radius  of  the  plasma  is  doubled  within  2pseconds. 

(4)  The  air  plasma  expands  from  a radius  of  . 2 cm  to  a radius  of  . 5 cm 
within  5 microseconds.  (5)  The  total  time  from  laser  turn  onto  production 
of  an  air  plasma  one  centimeter  in  diameter  is  therefore  of  the  order  of 
44  microseconds,  which  is  in  agreement  with  the  data  of  Ref.  4. 

Even  though  the  model  is  consistent  with  experiment,  there  remain 
several  internal  inconsistencies  which  should  be  eliminated  if  precise 
predictions  are  to  be  made.  For  example,  the  vapor  production  and  vapor 
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heating  are  clearly  coupled  aince  both  proceasea  operate  over  the  aame 
time  scales.  Similarly  the  final  atage  of  vapor  heating,  the  creation  of 
an  air  plaama  and  the  aubaequent  expanaion  of  the  air  plasma  all  proceed 
at  aimilar  ratea  and  ahould  be  united  into  one  continuous  process. 

Further  experimental  data  would  be  helpful  for  making  refinements 
to  the  model.  Experiments  ahould  be  carried  out  with  beams  which  are 
uniform  over  the  aiae  of  the  vapor  cloud,  say  . 2 cm.  If  the  intensity  does 
not  remain  constant  throughout  the  pulse,  the  important  quantity  for 
breakdown  calculations  is  the  fluence  at  breakdown,  not  the  time. 
Measurements  made  on  a variety  of  different  size  alumina  particles  having 
radii  greater  than  30p  would  aid  in  determining  whether  the  observed 
threshold  is  caused  by  jl  requirement  that  the  vapor  be  pure  or  a require- 
ment that  the  cloud  be  a certain  size.  The  threshold  intensity  should  not 
change  with  increasing  particle  size  if  the  threshold  is  related  to  the  purity 
of  the  cloud.  On  the  other  hand,  if  the  threshold  is  related  to  the  size  of 
the  cloud,  a larger  particle  will  have  a lower  threshold.  The  composition 
of  the  luminous  plasma  can  be  investigated  by  observing  atomic  line 
spectra.  The  fluence  at  which  the  luminous  vapor  first  appears  reveals 
information  about  the  particulate  vaporization  time.  Finally,  of  course, 
comprehensive  experiments  should  also  be  conducted  with  different 
particulates  and  at  different  wavelengths. 

In  conclusion,  we  have  shown  that  the  long  pulse  breakdown  of  air 
containing  30  micron  alumina  particles  can  be  explained  by  a model  based 
upon  laser  heating  of  the  particulate  vapor  with  the  subsequent  transfer  of 
energy  from  the  heated  vapor  plasma  to  the  surrounding  air.  However  the 
mechanism  depends  upon  several  details  of  the  laser  - particulate  interac- 
tion which  cannot  be  easily  extrapolated  to  different  wavelengths  or  to 
particulates  having  a different  composition. 
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